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FOREWORD 


The  Laser  Physics  Branch  of  the  Optical  Sciences  Division 
Naval  Research  Laboratory,  Washington,  D.C.,  prepared  this  Semiannual 
Report  on  work  sponsored  by  the  Advanced  Research  Projects  Agency 
ARPA  Order  20S2.  Co-authors  of  the  report  were  J  R.  Airey,  O.C.  Barr 
N.  DJeu,  J.L.  Emmett,  J.  P.olzrichter,  T.  Kan,  J.P.  Letellier,  J.  Marzolf, 
J.M.  McMahon,  S.  Searles,  J.  Stregack,  J.  Trenholme,  W.S.  Watt, 

W.  Whitney. 


ABSTRACT 


The  ARPA- NRL  high  energy  laser  program  is  concerned  with  the  develop¬ 
ment  of  laser  technology  in  four  program  areas;  Chemical  Lasers, 

Electric  Discharge  Lasers;  High  Power  Glass  Lasers  and  New  Laser 
Techniques.  This  report  summarizes  the  progress  made  in  those  areas 
during:  the  last  half  of  FY  72. 


PROBLEM  STATUS 

This  is  a  semiannual  technical  report;  work  is  continuing. 
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CHEMICAL  LASER  PROGRAM 


1.  DF-CCL  Supersonic  Transfer  Chemical  Laser  (TCL) 

The  concept  of  a  supersonic  DF-CO,  chemical  laser  capable  of  cw 
operation  without  a  vacu>s  pump  has  been  discussed  in  the  preceeding 
semiannual  technical  report. 

During  the  past  six  months  complete  installation  and  modification 
of  the  AVCO  Mark  III  GDL  has  been  accomplished  to  permit  DF-CO2 
chemical  laser  operation.  Figure  1  shows  the  completed  device  fully 
assembled  in  our  new  chemical  laser  facility.  Not  visible  in  the 
photograph  are  the  bottled  gas  tank  blocks  and  exhaust  gas  scrubber 
installation  located  outside  the  building. 

At  this  writing  the  cold  flow  tests  and  burner  testing  have  been 
completed.  We  are  now  making  measurements  of  exhaust  gas  composition 
to  establish  the  correct  flow  conditions  in  preparation  for  the  intro¬ 
duction  of  fluorine  and  deuterium  gases.  When  these  preliminary 
steps  are  completed  we  will  begin  a  series  of  laser  power  measurements 
with  a  multimode  mirror  configuration.  A  multiple  hole  output  mirror 
and  calorimeter  configuration  will  be  employed  in  this  series  of  tests. 
It  appears  reasonable  to  expect  laser  power  tests  to  begin  about 
October  1,  1972.  Dr.  Cool  is  returning  to  Cornell  University  at  the 
end  of  August,  but  will  continue  his  association  with  this  project  as 
a  consultant  to  NRL. 

While  the  laser  installation  has  occupied  most  of  our  time 
during  the  past  six  months,  attention  has  also  been  given  to  potential 
problems  associated  with  the  formation  of  carbonyl  fluoride  (COF2)  in 
the  present  device.  Our  measurements  of  the  rate  of  deactivation  of 
C02(00o 1)  by  COF2  and  the  measurement  of  absorption  of  CO2  laser 
radiation  by  COF2  performed  at  the  Air  Force  Rocket  Propulsion  Labora¬ 
tory  have  permitted  our  estimate  of  the  expected  effects  of  COF2  on 
laser  performance.  We  expect  no  serious  handicap  from  COF2  formation 
provided  the  partial  pressure  of  COF2  does  not  exceed  about  2  torr 
which  would  appear  to  be  an  unreasonably  large  amount  of  COF2  in  the 
present  device. 

2.  CO  Laser  Research 


Flame  Lasers 


We  have  experimentally  investigated  the  flame  for 

possible  optical  gain  in  the  product  CO  molecules.  Acetylene  was 
chosen  because  it  is  known  to  produce  vibrationally  highly  excited 
CO  molecules  in  reaction  with  0  atoms.  It  has  the  advantage  over  CS„ 
of  being  a  gas  at  room  temperature,  so  that  one  can  more  readily  go 
to  higher  pressures  with  the  ^2^2~^2 
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The  laser  oscillation  ringe  technique  was  used  to  probe  the 
C2H2-02  flame.  This  is  the  very  sensitive  method  for  measuring  gain 
developed  by  us  that  led  to  the  successful  operation  of  the  first 
flame  laser.  Quantitative  measurements  cf  the  optical  gain  for  the 
P9_g(12)  transition  were  made  across  the  profile  of  the  C-H  -0,  flame 
at  several  pressures  around  10  torr.  The  highest  gain  measured  was 
0.057®  cm"l.  The  results  of  the  work  lead  us  to  conclude  that  a  CO 
laser  based  on  that  system  is  feasible.  However,  because  of  the 
relatively  low  gain  attainable  in  the  flame,  such  a  laser  would 
probably  have  a  low  chemical  efficiency. 

CO  Transition  Probabilities 


Although  a  knowledge  of  the  transition  probability  is  central 
to  the  understanding  of  any  laser  system,  these  quantities  have  never 
been  measured  for  the  CO  laser  transitions  because  of  experimental 
difficulties.  We  have  succeeded  in  obtaining  the  CO  vibrational 
transition  probabilities  indirectly  by  measuring  the  ratios  of  transi¬ 
tion  probabilities  for  adjacent  bands  and  fitting  them  to  a  theoretical 
form. 


The  ratios  were  obtained  by  correlating  the  rotational  gain 
distributions  of  adjacent  bands  using  the  fact  that  the  distribution 
for  each  band  is  characterized  by  the  corresponding  inversion  ratio  and 
the  temperature  of  the  medium.  The  experimental  results  show  that  the 
square  of  the  electric  dipole  matrix  element  varies  linearly  with  the 
upper  state  vibrational  quantum  number  to  within  a  few  percent  for  the 
CO  molecule  through  V  =  12. 

The  particular  method  that  we  have  developed  for  this  experiment 
is  applicable  to  the  determination  of  relative  transition  probabilities 
for  cascading  bands  in  any  molecules.  A  paper  describing  the  techniques 
and  the  results  obtained  for  CO  has  been  submitted  to  the  Journal  of 
Chemical  Physics  and  is  given  as  Appendix  A  to  this  section. 


SHORT  PULSE  CO;  MOLECULAR  LASER 
!•  Short  Pulse  Oscillator 


n  «  Mode  -locking  and  pulse  shuttering  have  been  used  to  generate  a 
loCkPrC°MifUl8e  5r0m  a”  ato°8Pherlc  pressure  TEA  laser.  The  mode 
Internal  lolfof^h  a?°u8to-°Pt*c  effect  germanium  to  modulate  the 
of  ah  n88  the  laser  at  the  cavity  mode  separation.  Problems 
of  AR  coating  damage  to  the  modular  element  due  to  the  high  fluxes  In 
the  cavity  have  been  eliminated  by  the  use  of  a  Brewster- ancle  Ge 
modulator.  Alternative  techniques  to  achieve  mode- locking  of  the 
laser  utlilzlng  saturable  absorbers  and  a  liquid  CS,  acoufto-optlc 
ceii  nave  not  been  satisfactory  to  date.  The  saturlble  absorbers 

CS  acous^o-  M6  el1t?er  u°°  l0“  °r  t0°  hlgh  a  8aturation  flux  and  the 
f 2  "°P  ®  cell  exhibits  very  high  attenuation  to  acoustic-drive 

frequencies  in  the  region  of  interest  below  100  MHz.  The  pulse  shutter 

Une  tri8  ^  P°Ckela  cel’  ^ven  by  a  transmiLion 

of  >  1  if  £  8ur  trl8gered  8Park  8®P-  nanosecond  pulse 

of  >  1  mJ  is  limited  at  the  present  time  by  the  low  gain  of  the 

ren?6  h  i!**118,1  TEA  18881  °8cillator-  This  oscillator  will  be 
futire  rill  lnmeter4lon&  distributed  discharge  oscillator  in  the  near 

“S 

optics.  We  are  confident  5-10  nj  in  a  nanosecond  pulse  can  be^ener ated . 
Double  Raman  Conversion 


„.3+  Raman  scattering  experiments  in  H*  to  shift  the 

experiment  £  the  band  °f  C°2  have  heen^under taken.  These 

experiments  were  hampered  by  the  necessity  to  share  the  use  of  the  YAG 

short-pulse  oscillator  driver  for  the  glass  laser  system  T^ese 

experiment.  .Ill  be  continued  vhen  .  ..co„d  .eper.te  .hurt  put.” 

oscillator  system  no.  being  sssembled  Is  finished.  Although  51 

con.er.!on  of  1.06  u  to  1.9  „  hse  been  echleved,  the  seconfcon.er.lon 

scat  ter  ^served.  Since  the  Raman  gain  for  the  second 

scattering  process  is  lower  than  the  first  step,  competing  scattering 
process.,  could  be  depleting  the  1.9  micron  pu^.  These  ®u«tlo“f  * 
will  be  investigated  more  fully.  M 

3.  TEA  Pre-Amplif iera 

n..io«TW°  v<\iume  TEA  discharges  were  investigated  as  short 

pulse  preamplifier  elements.  One  utilizes  two  solid  Rogowski  electrodes 
8ur face  Preionization  created  by  UV  emission  from  an  auxiliary 
ischarge  while  the  other  utilizes  a  segmented  cathode  with  surface 
preionization  created  by  corona  discharges  between  metal  blades  and 
t.lass  rods.  The  solid  electrode  UV  initiated  devices  achieved  higher 

thekdfvic8'  b6t^eJ  PUl8e  tC  PUl8G  reliability  “d  were  more  durable  than 
?  1  ™AC°ntalning  dielectrics.  The  one  meter  long  2  cm  x  2  cm 
aperture  TEA  amplifier  has  displayed  a  short  pulse  gain  of  ~  20 
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indicating  a  gain  coefficient  of  3.27.  cm" ' .  Slightly  higher  gains  can 
be  achieved  with  additional  energy  storage.  A  second  identical 
amplifier  unit  is  now  being  prepared  to  be  followed  shortly  by  a 
third  3  cm  x  3  cm  amplifier.  The  final  6  cm  x  6  cm  preamplifier  will 
be  a  set  oi  Lumonix  9  a  x  9  cm  0.5  meter  TEA  lasers  now  under  procure¬ 
ment  and  due  to  arrive  at  NRL  in  November.  The  preamplifier  chain 
will  then  be  an  all  metal  system  without  fault-prone  dielectrics  and 
will  be  a  reliable  and  very  rugged  system. 

A.  10  Liter  e  beam  Amplifier 

Intensive  efforts  were  undertaken  in  this  last  reporting  period 
to  investigate  the  suitability  of  the  cold  cathode  high  curr  Vv,  density 
electron  gun  as  the  pre-ionizer  for  the  10  liter  e"beam  amplifier.  A 
joint  series  of  experiments  were  conducted  with  Maxwell  Laboratories 
and  the  Institute  of  Fluid  Mechanics  in  Marseille  to  study  the  laser 
properties  of  a  10  liter  discharge  stabilized  by  a  10  cm  x  100  cm  cold 
cathode  gun  de/ice.  The  experiments  were  successful  and  demonstrated 
that  high  current  density  beams  can  pui>p  CO2  in  a  short-pulse  energy- 
storage  mode.  A  peak  gain  of  A. 97.  cm"^  was  achieved  in  a  3 :  1  mixture 

of  He:N2:C02  for  a  energy  deposition  of  226  J/t.  This  figure 
represents  an  instantaneous  energy  storage  of  5-6  joules  of  optical, 
10.6  u  energy  which  can  be  extracted  in  a  short  pulse.  A  paper 
describing  this  experiment  is  included  in  the  latter  part  of  this 
report.  A  complete  10  cm  x  100  cm  cold  cathode  e"beam  system  is  now 
in  the  final  stages  of  procurement.  It  is  slated  to  be  delivered 
to  NRL  in  mid-December. 

5.  Saturable  Absorbers 


While  a  suitable  gaseous  or  gaseous-mixture  absorber  which 
suppresses  the  entire  CO2  gain  bands  has  not  yet  been  discovered, 
these  experiments  are  being  continued.  Present  studies  are  now  being 
concentrated  on  pressure  broadened  ^F,  and  mixtures  containing  this 
gas  because  it  displays  an  extremely  wide  absorption  band.  Since  a 
suitable  gas  mixture  which  can  cover  the  entire  CO2  spectrum  will  be 
difficult  to  find,  experiments  are  also  underway  to  investigate  a 
combination  of  a  dispersive  element  like  a  diffraction  grating  and  a 
passive  absorber.  The  diffraction  grating  would  select  out  a  single 
CO2  transition  by  angular  dispersion  along  the  optical  train  while  the 
saturable  absorber  would  be  needed  to  control  the  gain  in  the  single 
transition  of  interest. 
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NEW  LASER  TECHNIQUES 

During  this  reporting  period  several  experiments  have  been  planned 
and  their  construction  initiated. 

Prior  investigations  of  the  properties  of  potential  vapor  dye 
laser  molecules  were  conducted  using  a  laser  fluorescence  technique. 

A  N^  laser  was  borrowed  for  this  task  but  had  to  be  returned.  Conse¬ 
quently  the  purchase  of  a  N  laser  has  been  negotiated  and  it  will  be 
delivered  to  NRL  in  October  1972.  The  laser  will  be  used  to  help 
evaluate  lifetimes  of  potential  visible  laser  molecules. 

Several  optical  components  have  been  purchased  to  aid  in  the 
evaluation  of  solid  state  materials  for  high  energy  storage.  These 
include  etalons  and  frequency  doubling  crystals  for  YAG  systems.  This 
topic  has  a  direct  bearing  on  the  solid  state,  ARPA  supported,  program. 


Most  of  the  current  effort  has  been  directed  towards  the  production 
of  visible  laser  systems,  either  chemically  or  electrically  excited. 

It  is  unlikely  that  in  general  exothermic  chemical  reactions  will 
selectively  produce  electronic  state  inversions  in  product  molecules. 
Available  data  is  sparse  but  suggests  that  reactions  of  the  type 
M  +  x2  -  MX  +  X  hold  some  promise.  (M  -  rare  earth;  X  -  halogen). 
Fluorescence  spectroscopy  of  Zare  and  coworkers  on  the  reaction 
Ba  +  F2  "*  BaF  +  F  ha»  shown  that  the  A  and  B  states  of  BaF  are  not 
produced  but  that  the  higher  C  state  is.  This  is  perhaps  the  only 
example  of  a  simple  metathetlcal  reaction  where  an  electronic  inversion 
is  created.  In  order  to  capitalize  on  this  fact  a  reaction  system  has 
to  be  built  so  that  the  reagents  can  be  mixed  at  a  sufficiently  high 
density  to  produce  a  useful  gain  on  say  the  C  -  A,  or  C  -  X,  transitions. 
An  R-F  induction  heater  has  been  built  which  will  be  used  to  vaporize 
Ba  (or  Ca  etc.)  in  an  oven  which  will  act  as  a  high  pressure  source  of 
the  metal  atoms.  The  oven  is  also  presently  under  construction.  These 
will  form  a  part  of  a  high  pressure,  fast  flow,  fast  mixing  system  in 
which  the  alkaline  earth-halogen  reactions  will  be  investigated. 

_  ia  that  electrical  excitation  of  electronic  states  holds 

more  promise  in  the  future  than  chemically  excited  systems.  The 
electron  energy  in  a  system  can  be  tailored  to  match  specific 
excitation  cross-sections  and  the  range  of  energies  available  is  wider 
than  those  available  in  purely  chemical  driven  schemes.  The  concept  of 
double-pulse  electrical  excitation  is  being  pursued.  For  example,  all 
attempts  to  initiate  copper  vapor  laser  action  using  a  simple  pulse 
technique  using  a  copper  compound  as  the  starting  material  rather  than 
atomic  copper  have  30  far  failed.  The  discharge  conditions  in  a  single 
pulse  are  most  probably  not  suitable  for  both  dissociation  and  subsequent 
excitation.  A  double-pulse  technique  where  each  pulse  may  be  engineered 
to  perform  the  individual  tasks  should  hold  more  promise.  Similarly,  if 
laser  action  is  to  be  attempted  using  chemically  metastable  radicals  as 
the  medium,  it  would  be  useful  to  have  a  pulse  to  form  the  radicals 
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(which  takes  a  finite  time  governed  by  the  reaction  kinetics) 
by  the  ejccitation  pulse.  A  combined  R-F,  TEA  double  discharge 
apparatus  is  being  designed  to  use  as  indicated  { bove  8 


followed 

pulse 
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HIGH  POWER  GLASS  LASER  PROGRAM 


1.  Introduction 


At  the  close  of  the  last  reporting  period  there  were  clearly  a 
number  of  problem  areas  which  required  resolution  if  the  NRL  laser 
system  was  to  operate  at  or  near  specification: 


(a)  the  original  configuration  of  the  French  VD  640  laser 
appeared  drastically  inappropriate  for  short  pulse  amplification, 
ntensive  research  effort  was  launched  to  optimize  the  system  for 
subnanosecond  pulses; 


An 


(b)  a  parasitic  oscillation  problem  was  discovered  in  the  disc 
amplifier  caused  by  self  oscillation  of  individual  discs  at  stored 

*evela  °f  about  °’4  J/«n3.  An  experimental  program  was  begun 
in  which  we  evaluated  a  series  of  oscillation  suppressing  edge 
coatings  developed  by  Owens  Illinois  Glass  Co. 

A  parallel  but  closely  coupled  design  effort  to  design  a  1-2  kJ 
hand  amplifier  has  been  Punched  and  a  preliminary  design  is  in 


We  have  also  developed  a  mode- locked  Nd:LaCaSoAp  oscillator 
which  produces  time-bandwidth- limited  operation  for  pulses  as  short 
as  5  psec.  This  represents  a  substantial  extension  of  the  stat\*  of 
the  art  from  mode  locked  Nd: glass  oscillators. 


2.  French  Laser  Modifications 

During  operation  of  the  VD  640  with  subnanosecond  pulses  in  the 
previous  reporting  period,  self- focusing  damage  was  noted  at  levels 
well  below  that  expected  from  reports  of  short  pulse  operation  at 
other  laboratories  (Hindsight  reveals  that  in  the  cases  of  Sandia 
and  Rochester  a  large  fraction  of  the  energy  was  not  contained  in 
the  short  pulse  but  In  the  background.)  Consideration  of  the  self- 
jcusing  problem  revealed  several  generally  true  prescriptions  for 
operation  when  self- focusing  is  the  limiting  damage  mechanism: 

-  use  a  minimum  path  length  of  glass  (i.e.  maximum  gain  per 
unit  length), 

-  propagate  a  divergent  spherical  wave  through  the  system  as 
this  will  have  a  much  higher  self- focusing  threshold, 

-  propagate  a  TEM  spatial  mode. 

oo 

Experiments  on  self- focusing  in  the  various  laser  glasses  were 
performed  and  it  was  found  that  ED-2,  LG-56  and  MG-915  had  essentially 
the  same  nonlinear  index  of  refraction  and  additionally  the  non¬ 
linear  index  of  refraction  was  a  constant  for  pulses  from  20  psec 
to  250  psec.  This  implies  that  the  self- focusing  can  follow  any 


8 


temporal  modulation  on  the  pulse  and  give  a  much  lower  threshold  for 
self-focusing.  Thus,  a  further  prescription  for  optimum  operation  is 

-  use  a  time-bandwidth  limited  oscillator. 

These  ^periments  are  reported  in  more  detail  in  Appendix  1. 

These  dicta  were  implemented  in  the  following  fashion  to  optimize 
the  VD  640. 

-  modification  of  the  oscillator  to  incorporate  a  bandwidth 
selecting  etalon  as  the  front  mirror  of  the  oscillator.  In 
this  mode  of  operation,  if  a  pulse  is  generated  it  has  no 
fine  grain  modulation  (if  the  etalon  is  misaligned  no  pulse 
is  generated) .  Zero  degree  sapphire  was  chosen  since  the 
etalon  reflectivity  (25%)  is  very  near  the  optimum  coupling 
reflectivity, 

-  Owens-Illinois  ED-2  rods  were  obtained  for  the  system  (on 
another  contract).  The  higher  gain/ joule  stored  of  this 
glass  allowed  us  to  reduce  the  path  length  of  glass  in  the 
system  by  85  cm  and  also  to  stage  for  a  4  x  energy  gain/stage, 

-  very  careful  alignment  procedures  were  worked  out  to  ensure 

that  the  mode  propagated  would  in  fact  remain  very  close  to 

a  TEM  spatial  mode, 
oo 

Figure  1  shows  the  present  system  configuration  as  used  for 
experiments  (through  the  45  mm  stage)  and  for  disc  laser  testing 
with  the  full  system.  The  improved  alignment  has  worked  quite  well 
to  a  point,  in  that  through  the  45  mm  amplifier  25  joule,  250  psec 
and  40  joule,  900  psec  pulses  have  been  obtained  with  an  essentially 
diffraction  limited  beam  quality  for  several  hundred  shots.  Figure  2 
shows  a  shear  plate  interferogram  of  the  pulse  at  25  J  in  250  psec 
(10  watts)  before  recollimation.  It  can  be  seen  that  the  maximum 
deviation  from  fringe  linearity  is  ~  1/3  fringe. 

Unfortunately,  this  happy  state  of  affairs  does  not  prevail 
through  the  64  mm  rod  at  full  power  when  there  is  pre-existing 
internal  (or  surface)  rod  damage.  Up  to  =-  40-50  joules  in  250  psec 
the  beam  is  well  behaved  and  the  shear  pattern  shown  in  Figure  3a 
is  obtained.  Above  that  level  a  rather  dramatic  change  takes  place. 

The  beam  essentially  blows  up  in  the  sense  that  the  divergence 
becomes  abnormally  high  (~  20-30  mr)  and  the  beam  becomes  essentially 
depolarized  and  incoherent.  Figure  3b  shows  a  shear  pattern  taken 
when  this  effect  is  occurring.  Experiments  were  run  to  elucidate 
the  functional  dependence  of  this  phenomena  at  250  psec  and  900  psec. 
These  experiments  substantiated  several  features: 

-  when  the  effect  occurs  the  output  becomes  essentially  incoherent 

in  that  when  it  intercepts  a  straight  edge,  no  diffraction 

pattern  downstream  is  created. 


S 


the  effect  seems  to  scale  with  power  density  in  that  the  onset 
is  proportionally  higher  at  longer  pulsewidths  and  near 
threshold  only  the  beam  center  will  show  this  behavior, 

geometrical  effects  (thermal  lensing  in  the  laser  rods)  are 
not  of  primary  importance, 

the  effect  is  not  gross  self- focusing  of  the  entire  beam  in 
that  the  output  wave  does  not  originate  from  a  point,  but 
rather  corresponds  to  severe  wavefront  distortion  on  a  point- 
by-point  basis, 

the  effect  is  related  to  the  past  history  of  the  laser  rod  in 
that  a  rod  with  no  damage  will  operate  without  showing  this 
effect  up  to  quite  high  levels  (>  100  J  in  250  psec). 
Realistically  it  is  a  limitation  to  routine  operation,  since 
damage  will  occur  sooner  or  later. 


The  most  likely  candidate  effect  is  wavefront  distortion  around 
any  damaged  regions  in  the  rod  which  at  high  intensity  is  strong 
enough  to  stabilize  the  Fresnel  diffraction  pattern  downstream  from 
the  point  of  inception  and  give  very  severe  localized  wavefront 
distortion  and  subsequent  beam  steering.  This  mechanism  has  the 
requisite  features  to  explain  the  experimental  observations  and  the 
level  at  which  it  occurs  (~  1010  watts/ cm2)  is  reasonable. 

3.  Disc  Laser  Test  Results 


Since  the  close  of  the  last  reporting  period,  the  parasitic 
osciliation  problem  in  the  present  disc  laser  has  been  solved  and  a 
sufficient  theoretical  understanding  has  been  developed  to  define  the 
requirements  for  edge  coatings  for  larger  disc  lasers. 


To  achieve  the  first  objective,  a  black  coating  was  designed  by 
Owens-Illinois  and  applied  to  the  NRL  laser  discs.  The  coating  had 
an  index  of  refraction  of  =“1.8  at  1.06  u  and  was  sufficient  to 
suppress  parasitic  oscillation  at  energy  storages  up  to  the  design 
level  of  0  5  I/cm-5  (see  Appendix  B  for  coating  test).  A  series  of 
experiments  on  the  small  signal  gain  coefficient  of  the  disc  laser 
with  coated  and  uncoated  discs  were  performed  (Appendix  C)  and  the 
performance  with  the  coated  discs  was  found  to  agree  closely  with  the 
computer  calculations. 


The  present  generation  coating  is  not  unambiguously  a  success 
in  that  the  laser-glass  coating  interface  must  be  shielded  from  direct 
lamp  irradiation  to  prevent  shear  damage  In  the  laser  glass  caused 
y  differential  heating.  Shields  have  been  designed  and  installed 
which  prevent  this. 


uTaj}6  I, lists  the  results  the  latest  sequence  of  active  tests 
of  the  disc  laser.  It  was  not  possible  to  achieve  operation  in  the 
400-500  joule  regime  because  of  the  effects  noted  in  the  previous 
section  but  an  output  of  320  joules  in  a  900psec  pulse  was  obtained. 
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The  chief  hinderance  here  was  some  depolarization  of  the  input  beam 
leading  to  a  loss^  coefficient  of  =-  27.  cm,  rather  than  the  much  lower 
value  of  .5  -  .75%/cm  obtained  with  a  well  polarized  input.  The  table 
lists  the  pump  clergy,  input  energy,  actual  output,  computed  output 
(with  the  actual  loss)  and  computed  output  with  the  lower  loss. 


TABLE  I 


Ep(kJ) 

ein<j> 

eout^U; 

Wn> 

ecalc^*75^ 

118 

60 

220 

220 

250 

113 

70 

280 

275 

320 

150 

110 

320 

325 

395 

The  disc  edge  shields  occluded  207.  of  the  disc  area  so  that  a 
smaller  beam  diameter  than  that  originally  envisioned  had  to  be  used 
(50  cm  rather  than  60  cm2)  which  lowered  the  output  level  by  =«  10%. 

The  output  beam  appeared  to  be  well  behaved  and  not  subject  to  the 
beam  breakup  phenomena  noted  with  the  final  rod  amplifier  although 
the  energy  density  was  50%  higher.  This  is  reasonable  since  the  path 
length  through  a  disc  is  very  short  compared  to  an  equivalent  rod 
amplifier. 

In  summary,  the  disc  amplifier  seems  to  functionally  behave  as 
designed  and  we  feel  that  the  physics  questions  as  to  the  viability 
of  this  concept  have  been  answered.  We  are  still  limited  at  the 
present  time  ^n  preparing  a  suitable  input  pulse  to  run  the  system 
to  energies  in  exces  of  40G  joules  but  at  least  have  a  good  under¬ 
standing,  of  system  behavior  at  energy  densities  equivalent  to 
operation  at  400  joules. 

4.  1-2  kJ  Disc  Amplifier  Design 

During  the  reporting  period  several  design  alternatives  for  a 
1-2  kJ  disc  amplifier  were  generated  and  tradeoffs  done  on  each  one. 
The  design  constraints  were: 

-  output  energy  of  at  least  1  kJ  with  a  300  joule  input  pulse; 
output  of  at  least  2  kJ  with  a  500  joule  input  pulse, 

output  energy  density  no  higher  than  in  the  present  device 
at  the  500  joule  level, 

-  aperture  size  no  larger  than  practical  to  fabricate  from  a 
glass  manufacturing  standpoint  at  the  present  time, 

-  gain  across  the  aperture  no  larger  than  practical  from  a 
parasitic  oscillation  standpoint  (see  Appendix  D) , 
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-  coupling  efficiency  assumed  to  follow  ZAP  scaling  as  lamp 
packing  densities  are  changed, 

-  loss  coefficient  assumed  to  be  ~  1%/cm. 

Three  design  alternatives  were  chosen  and  the  major  features  of 
each  are  specified  in  Table  II.  The  first  entry  is  a  direct  scaleup 
of  the  present  system  to  twice  the  aperture  size  and  the  last  repre¬ 
sents  the  largest  disc  Owens-Illinois  feels  it  is  practical  to 
fabricate  with  good  optical  quality  at  this  time.  The  intermediate 
case  represents  a  system  with  an  output  intensity  at  2  kJ  equivalent 
to  what  we  have  run  in  the  present  disc.  In  all  cases  the  energy 
storage  was  varied  to  give  a  gain  no  higher  than  exp.  (2.5)  across 
the  disc. 

Present  indications  are  that  either  the  first  or  the  second  (and 
probably  the  second)  is  the  best  choice  from  the  standpoint  of 
efficiency.  One  further  design  constraint  has  been  added  from  a 
systems  standpoint: 

-  the  aperture  size  shall  be  small  enough  to  permit  a  Faraday 
rotator  isolator  to  be  added  on  the  output  to  isolate  the  laser 
from  the  target. 

The  tradeoff  here  is  not  as  clear  in  detail  but  certain  features 
are  obvious.  For  a  given  uniformity  (i.e.  extinction  ratio)  the 
energy  in  the  rotator  magnetic  field  will  increase  as  the  cube  of 
the  diameter.  The  difficulties  of  obtaining  high  Verdet  constant, 
high  damage  threshold  glass  in  this  size  cannot  be  exactly  evaluated 
yet  but  they  are  appreciable.  Table  III  illustrates  this  dilemma 
for  a  coil  designed  for  1GJ  extinction  ratio  with  Owens-Illinois 
EY-1  and  Schott  SF-6 

TABLE  III 

DESIgN  »Wr,kJ  e(sf-mu 


1 

60 

360 

2 

130 

780 

3 

200 

1200 

Owens-Illinois  EY-1  is  a  variant  of  the  basic  ED-4  glass  with 
40%  Terbium  oxide  by  weight.  If  the  glass  technology  proves  out  this 
glass  will  be  the  choice  for  large  aperture  Faraday  rotators.  The 
numbers  given  in  Table  III  are  based  on  the  measurements  done  at  NRL 
(Appendix  E). 


13 


5.  Short  Pulse  Bandwidth  Limited  Oscillators 


During  this  reporting  period,  we  have  investigated  extending  the 
methodology  used  to  develop  the  mode  locked  Nd:YA<tG  lasers  to  other 
hosts  with  the  objective  of  generating  shorter  pulses  while  maintaining 
time-bandwidth  limited  operation. 

The  selection  criteria  for  materials  for  this  usage  *?ere: 

-  broader  fluorescence  linewidth  than  neodymium: YAG 

-  crystalline  host  so  that  the  nonlinear  index  of  refraction 
would  be  small  compared  to  Nd: glass. 

The  first  candidate  material  tried  has  been  evaluated  and  found 
to  perform  with  fairly  spectacular  results.  This  material  is  Nd  doped 
into  a  Lanthanum  Calcium  Sodium  Apatite  structure  (SOAP)  which  has 
been  developed  by  Westinghouse.  Our  results  to  date  indicate:  (see 
Appendix  F) 

-  time  bandwidth  limited  operation  ~  5-6  psec  pulses 

-  pulse  trains  with  >  107.  of  the  energy  in  the  strongest  pulse 

-  an  emission  wavelength  10611. 4A  which  overlaps  the  second 
strongest  transition  in  Nd:YA-tG  preamps  which  can  be  used  to 
amplify  the  oscillator  pulses.  A  formal  publication  is  in 
preparation  at  this  time.  Figure  4A  shows  a  typical  pulsetrain 
and  Figure  4B  shows  the  output  spectrum  of  the  pulses. 


aser  Configuration 
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APPENDIX  a 

A  METHOD  FOR  MEASURING  RELATIVE  TRANSITION  PROBABILITIES 
OF  CASCADING  MOLECULAR  BANDS: 

APPLICATION  TO  CO  FUNDAMENTAL  BANDS’’ 

N.  Djeu 
S.K.  Searles 

Laser  Physics  Branch 
Optical  Sciences  Division 

Naval  Research  Laboratory 
Washington,  D.C.  20390 

ABSTRACT 


A  general  method  is  presented  to  determine  the  relative  transition 
probabilities  for  two  molecular  bands  which  have  one  vibrational 
level  in  common.  It  is  based  on  optical  gain  (or  loss)  measurements  on 
a  number  of  vibration-rotation  transitions  in  each  of  the  two  bands. 

The  method  has  been  applied  to  CO  fundamental  bands  in  the  ground 
electronic  state  to  yield  |R^  |^/|r^_^|2  for  v  *  6  through  v  ■  11, 
where  Rv  ^  is  the  rotationless  matrix  element  of  the  dipole  moment 
connecting  the  v-lth  and  vth  vibrational  levels.  The  results  were 
fitted  by  using  the  wave  function  approximation  for  R^_ .,  and  the  ratio 
Rp  /Rq1  -  0*0002^^  was  obtained. 


I.  INTRODUCTION 

Experimental  radiative  transition  probabilities  are  directly 
available  only  for  the  lowest  few  vibrational  bands  in  molecules,  since 
they  are  determined  from  absorption  or  emission  intensity  measurements 
which  require  appreciable  known  populations  in  the  absorbing  or 
emitting  states.  Relative  transition  probabilities  can  be  found  for 
bands  which  share  the  same  upper  or  lower  vibrational  level,  inasmuch 
as  no  knowledge  of  the  population  density  of  the  common  level  is 
required.  But  conventional  intensity  measurements  do  not  yield 
information  on  the  relative  strengths  of  cascading  bands.  In  this 
paper  we  report  a  general  method  which  gives  the  transition  probability 
for  one  band  relative  to  that  of  a  second  which  originates  from  the 
vibrational  level  that  the  first  band  terminates  in.  This  opens  up 
the  possibility  of  determining  the  absolute  transition  probabilities  for 
highly  excited  bands  from  known  low  lying  band  intensities  through 


*  This  work  was  supported  under  ARPA  Order  2062  by  the  Advanced 
Research  Projects  Agency. 
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successive  ratios. 


In  the  case  of  CO,  the  v  *  1  ■*  v  ■  0  in  the  ground  electronic 
state  is  the  only  fundamental  band  for  which  the  transition  probability 
has  been  measured^.  Though  various  calculations  yield  values  for  the 
squares  of  the  electric  dipole  matrix  elements  which  agree 

with  each  other  to  within  a  few  percent  up  through  v  10,  there  has 
never  been  any  direct  experimental  evidence  to  support  those  calcula¬ 
tions.  In  measuring  directly  the  ratios  |2/  |R^  P,  the  present 

experiment  corroborates  the  results  of  those  calculations.  The  ratios 
i»r1|2/K.,i2  are  first  obtained  from  gain  measurements  on  the  P 
branch  transitions  of  the  v  *  6  H  v  ■  5  through  v  ■  12  ”•  v  *  11  bands 
in  a  cw  CO  laser.  They  are  then  fitted  to  the  theoretical  model  using 
the  wave  function  approximation  for  the  matrix  elements.  Each  R^_.  1® 
expanded  in  terms  of  Ro*,  with  the  expansion  truncated  after  the  first 
four  terms.  Using  known  values  of  Rq  /Rq  and  Rq3/Rq*,  We  determine 
RoVRq1  from  the  experimental  data.  UThese  numbers  are  then  used  to 
calculate  |R^_^ |2  in  terms  of  (Rq1!2. 

Because  of  the  relatively  low  gain  present  in  the  CO  laser,  we  had 
to  resort  to  a  rather  unconventional  method  for  gain  measurements  in 
order  to  obtain  highly  accurate  data.  A  single  mode,  frequency  scanned 
laser  is  used  for  that  purpose.  The  laser  oscillation  range  technique 
for  measuring  relative  optical  gains  is  described  in  some  detail  in 
Section  III.  Additional  considerations  about  this  technique  can  be 
found  in  a  forthcoming  paper3. 

II.  BASIS  FOR  EXPERIMENT 

Consider  a  molecular  system  with  its  vibrational  levels  greatly 
overpopulated  with  respect  to  the  translation-rotation  temperature. 

These  conditions  are  commonly  found  in  electrical  discharges,  exothermic 
chemical  reactions,  gasdynamic  expansions,  etc.  The  optical  gain  (or 
loss)  distribution  of  the  rotational  fine  structure  of  a  vibrational 
band  in  such  a  system  is  characterized  by  the  rotational  temperature 
and  the  ratio  of  the  population  densities  of  the  two  connected 
vibrational  levels.  The  distribution  is  a  very  sensitive  function  of 
the  latter  if  the  two  vibrational  levels  have  comparable  population 
densities.  Thi H  fact  allows  one  to  determine  the  relative  densities 
of  the  two  vibrational  levels  from  the  rotational  gain-loss  distribution 
of  the  band.  If  measurements  are  made  on  two  cascading  bands  in  a 
molecule,  the  ratio  of  the  two  corresponding  transition  probabilities 
can  be  found.  In  what  follows  we  restrict  our  attention  to  diatomic 
molecules,  but  the  principle  can  clearly  be  extended  to  polyatomic 
molecules. 

The  small  signal  optical  gain  coefficient  for  the  P  branch 
vibration-rotation  transition  (v,  J-l)  —  (v-1,  J)  with  mixed  Doppler 
and  pressure  broadening  can  be  written  as  follows:  * 
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K-ii2  J  /<v>p(j)> 


v,v-i,  % 

°p(j)  <«> 


["vBv  «p(-rv<J-»h  <=/kT)  -  »v.lVl  «p(-Fv.l<J)h  c/kt)] 

x  Re  v(2i*#lfW)'Fi  + 1  ^K)J>)  Vs"1) 


a) 


Here  k  Is  the  Boltzmann  constant,  T  is  the  temperature  of  the  medium 
(assumed  to  be  the  same  for  translation  and  rotation),  M  is  the  molecular 
mass,  c  is  the  velocity  of  light,  and  N  and  B  are  the  population 
density  and  the  rotational  constant  forVthe  vtK  vibrational  level.  Also 
lRy_ll  and  P(J))  are  the  rotationless  matrix  element  and  the 

vibration-rotation  interaction  factor  for  the  transition,  and  v(v,  P(J)) 
and  tv (v,  P(J))  are  the  line  center  frequency  and  the  pressure  broad¬ 
ened  width  (including  natural  broadening)  of  the  transition.  The 
rotational  energy  function  F  (J-l)  and  the  Doppler  broadened  width 
AVp(v)  are  given  by  v 


Fv(J-l)  -  Bb(J-l)  J 

and 

h(v)  =  t*  * 2- 


Finally,  Re  w  (x  +  iy)  is  the  real  part  of  the  complex  error  function. ^ 

This  somewhat  involved  expression  can  be  understood  in  two  parts. 
The  part  preceeding  Rg  w  gives  the  gain  at  linecenter  of  the  vibration- 
rotation  transitions,  while  the  function  Re  w  itself  is  simply  the 
lineshape  function  for  mixed  pressure  and  Doppler  broadening.  Note 
that  the  linecenter  gain  distribution  is  not  appreciably  altered  by 
pressure  broadening  so  long  as  the  J  dependence  of  the  pressure 
broadened  width  is  slight. 

In  general,  all  the  spectroscopic  constants  in  Eq.  (1)  are  known 
so  that  the  only  parameters  to  be  determined  are  T,  tv  ,  |RV  |2  N  and 
Vi-  However  the  last  three  quantities  only  appear  ifi  the” form  of  the 
products  |Ry_^l  and  |R  , |  N  . .  This  means  that  there  are  in 
effect  only  four  unknowns  in  Eq.  (I).  It  also  means  that  neither  the 
matrix  element  jR^  |2  nor  the  Population  densities  Nv  and  Nv  l  can  be 


A-3 


determined  uniquely  from  fitting  Eq.  (1)  to  experimental  data  on 
rotational  gain  distributions.  However,  if  gain  distributions  are 
known  for  two  successive  bands  v  +  1  -•  v  and  v-*v  -  1,  it  is  clear 
that  tjje  ratio  |r™ 12/  |Rj_ ,  |2  can  be  obtained  from  Jr^1  |2  N  and 
'Rv-1'  Ny‘  Furthermore,  only  relative  gains  for  the  two  bands  are 
needed,  since  the  population  density  of  the  common  vibrational  levels 
cancels  itself  when  the  ratio  is  taken. 

For  some  applications,  the  quantity  of  real  physical  significance 
is  the  Einstein  A  coefficient.  For  the  PfJ)  transition  it  is  related 
to  the  rotationless  matrix  element  |r^  |2  by  the  equation 


4  K-i!2  J  /(v>  P(J» 


A(v,  J-l)  - 


64  it 
3h\J 


2  J-l 


(2) 


So  long  as  the  vibration-rotation  interaction  is  small  (i.e.  /(v,  P(J)) 
=•  1)  ,  we  have 


ACvfl.  J-l) 
A(v,  J-l) 


(im\3  I Cf 
\W  K-;!2 


(3) 


v 

where  Xp.j.  is  the  wavelength  of  the  P(J)  transition  in  the  v  -  v  -  1 
band.  It  Should  be  pointed  out  that  the  ratio  in  (3)  is  taken  between 
A  coefficients  for  transitions  connecting  the  same  rotational  levels. 
This  condition  is  implicitly  assumed  whenever  we  mention  ratio  of 
Einstein  A  coefficients  in  this  paper. 

III.  TECHNIQUE 

The  ratios  |R  |  /  |r  ,  |  are  expected  to  vary  only  slightly 
with  v  after  the  first  few  vibrational  levels.  It  is  therefore 
essential  that  highly  accurate  relative  gain  distributions  be  obtained, 
so  that  the  v  dependence  of  the  matrix  element  ratios  is  not  obscured 
by  the  scatter  of  the  experimental  results.  Conventionally  gain 
measurements  are  made  by  sending  the  output  from  a  laser  oscillator 
through  a  second  gain  tube  (the  amplifier)  and  noting  the  increase  in 
the  transmitted  signal  when  the  amplifier  is  turned  on.  But  with  that 
method  it  is  very  difficult  to  measure  small  pains  with  good  accuracy. 
Typically  the  relative  error  for  a  measurement  giving  a  107.  optical 
gain  would  be  roughly  107..  To  circumvent  this  difficulty,  we  have 
developed  a  new  technique  which  has  the  additional  advantage  of 
requiring  only  one  gain  tube,  the  oscillator  itself.  The  remainder  of 
this  section  will  be  devoted  to  a  detailed  description  of  this  new 
technique . 


We  begin  with  a  discussion  on  the  operation  of  a  single  mode 
frequency  tuned  laser.  Consider  first  the  single  mode  oscillation 
condition  in  a  laser  illustrated  in  Fig.  1.  The  gain  profile  of  the 
laser  medium  as  a  function  of  frequency  is  given  by  the  factor  Re  w  in 
Eq.  (1).  The  cavity  loss,  on  the  other  hand,  is  generally  frequency 
independent.  It  is  represented  by  the  horizontal  dashed  line.  The 
laser  cavity  as  a  whole  has  net  gain  only  in  the  frequency  range  Av. 

But  the  condition  that  there  is  net  gain  over  part  of  the  gain  curve 
does  not  by  itself  mean  that  there  will  be  laser  oscillation.  Whatever 
medium  is  inside  a  laser  cavity  of  length  £,  there  can  be  oscillation 
only  at  the  discrete  set  of  frequencies  given  by  n  c/2£,  where  n  is  an 
integer. “  The  cavity  is  said  to  have  modes  separated  by  c/2£  shown 
as  short  vertical  lines  on  the  frequency  axis.  Whenever  one  of  these 
cavity  modes  falls  at  a  frequency  for  which  laser  gain  exceeds  cavity 
loss,  there  will  be  oscillation  at  that  frequency.  For  the  case 
illustrated  in  Fig.  1,  there  would  be  oscillation  at  the  frequency  of 
the  middle  mode.  Furthermore,  the  middle  mode  is  the  only  mode  that 
falls  within  Av,  creating  a  situation  known  as  single  mode  operation. 

When  the  cavity  length  is  changed,  each  cavity  mode  will  move 
along  the  frequency  axis.  In  that  process,  the  frequency  of  laser 
oscillation  changes  accordingly,  resulting  in  a  frequency  scanned  (or 
tuned)  laser  output.  For  a  cavity  length  change  equal  to  one  half  the 
wavelength  of  the  laser  transition,  one  cavity  mode  is  moved  to  the 
position  previously  occupied  by  its  adjacent  mode.  In  practice,  the 
change  in  cavity  length  can  be  accomplished  by,  for  example,  the  use 
of  a  piezoelectric  transducer  which  converts  an  applied  ramp  voltage 
into  a  linear  translation.  So  long  as  Av  <  c/2£  there  can  be  at  most 
one  mode  oscillating  at  any  time.  The  output  of  such  a  laser  is  then 
called  a  single  mode  frequency  tuning  curve.  Since  Av  is  the  frequency 
range  over  which  there  is  oscillation  as  the  cavity  length  is  swept, 
it  is  known  as  the  laser  oscillation  range.  A  schematic  of  a  single 
mode  tuning  curve  is  illustrated  in  Fig.  3  near  the  detector.  The 
frequency  range  covered  by  each  mode  is  c/2£  centered  at  the  line  center 
of  the  transition. 

Consider  now  a  distribution  of  gain  profiles  for  a  number  of  P 
branch  transitions  in  a  vibration-rotation  band  as  illustrated  in 
Fig.  2.  The  actual  spacing  between  lines  is  typically  10^  times  the 
linewldth  at  a  pressure  of  a  few  torr.  The  precise  distribution  of 
these  lines,  including  the  frequency  dependence  of  each,  is  given  by 
Eq.  (1).  In  the  particular  example  shown  in  Fig.  2  three  of  the 
transitions  have  sufficient  gain  to  overcome  threshold.  Therefore, 
each  of  these  can  be  made  to  produce  a  single  mode  frequency  tuning 
curve  in  a  wavelength  selective  cavity,  provided  that  the  oscillation 
range  for  the  P(J)  transition  AVp/.x  <  c/2£  for  each  J.  Thr.  dispersive 
element  used  for  wavelength  selection  is  typically  a  prism  or  a 
diffraction  grating.  If  the  wavelength  dependence  of  its  optical  loss 
is  slight  in  the  region  of  interest,  one  can  assume  that  the  cavity  loss 
is  the  same  for  all  transitions.  It  is  clear  then  that  there  is  a 


definite  relationship  between  the  cavity  loss  and  the  oscillation 
ranges  ^(j)  on  the  one  hand  and  the  parameters  in  Eq.  (1)  on  the 
other  hand.  Given  a  sufficient  number  of  Avp,m  for  a  fixed  cavity 
loss,  one  should  be  able  to  determine  the  products  |RV  . I2  N  and 

IJv-l'.  Nv-1  <relative  to  the  cavity  loss),  the  pressure  broadened 
wi ith  Av  ,  as  well  as  the  temperature  of  the  medium  T.  Suppose  one 
f  fra<luency  independent  variable  loss  l,  e.g.,  variable 
diameter  iris,  into  the  laser  cavity  in  addition  to  the  fixed  loss 
“*  Then  the  experiment  can  be  done  for  a  number  of  L  +  l  and 
statistical  averages  can  be  taken.  * 


IV.  EXPERIMENT  AND  RESULTS 

The  experimental  arrangement  used  for  making  the  relative  gain 
measurements  is  shown  schematically  in  Fig.  3.  The  CO  gain  tube 
contained  a  flowing  mixture  of  0.3  torr  CO,  0.2  torr  air,  and  2.5  torr 
?*  1  excited  by  a  3  mA  D.C.  glow  discharge  with  less  than  1% 

r  pple.  The  discharge  was  0.5  m  long,  and  approximately  85%  of  its 
length  was  cooled  by  liquid  nitrogen.  The  2.55  m  cavity  was  formed  by 
a  3  m  radius  of  curvature  silver  mirror  and  a  300  lines/mm  grating 
blazed  for  4  /an.  The  mirror-grating  assembly  was  mounted  on  an  air- 
sags  suspended  steel  beam  which  was  physically  isolated  from  the 
discharge  tube.  The  mirror  was  attached  to  a  piezoelectric  transducer 
modulated  by  a  ramp  voltage  repeated  at  ~  1  sec  intervals.  A  variable- 
iameter  iris  was  placed  directly  in  front  of  the  mirror  both  to  insure 
fundamental  mode  operation  and  to  act  as  a  variable  loss  device.  To 
m  nimize  water  vapor  interference,  the  regions  between  the  NaCl 
Brewster  s  angle  windows  and  the  mirror  (grating)  were  flushed  with 
N2  except  for  about  10  cm  near  the  grating.  The  frequency  scanned 
laser  output  was  coupled  off  the  grating  and  sent  into  a  gold  doped 
germanium  detector.  One  sweep  of  the  piezoelectric  transducer  would 
generally  cover  two  modes.  Photographs  were  taken  of  the  frequency 
scanned  signal  displayed  on  an  oscilloscope.  Since  the  separation 
between  adjacent  modes,  c/2£,  is  known,  the  oscillation  ranges  were 
readily  calculated. 

of  Wa?  d°ne  ln  tW0  8tePs‘  The  first  was  a  determination 

of  the  lineshape  function  Re  w,  for  which  Av  and  T  were  needed.  The 
oscillation  ranges  of  four  or  more  P  branch  transitions  in  one  band  were 
measured  for  a  given  setting  of  the  iris.  For  a  chosen  set  of  and 
T  the  line  center  gain  relative  to  total  cavity  loss  L  +  l  £as 

calculated  for  eac l  transition  according  to 
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where  the  J  dependence  of  Av^  has  been  dropped.  Then  a  linear  least 

square  fit  by  varying  |RV  ^  N  and  |rV  . |2  N  .  was  performed  to 
minimize  V  V  V  V 


exp 

°P(J)  “  ^CJ) 


exp 

^(J) 


(5) 


where  cl,, n  is  the  line  center  gain  calculated  from  Eq.  (1)  using  the 
assume?!,  |RV  .  |2  N  ,  and  |RV  ,  |2  N  _7.  The  value  of  v*.  minimized 

by  the  parameters  |r^  ^|2  and  |r^  ^ I  Nv_^,  as  a  function  of  Av^ 

and  T  showed  a  local  minimum  at  T  ■  90°K  and  Av  / Av^  *  0.27.^  These 
numbers  are  very  reasonable  in  view  of  the  following  considerations. 

The  value  of  the  optical  collision  frequency  of  6.2  x  10?  sec-1  implied 
by  the  measured  AVp  is  to  be  compared  with  the  hard  sphere  collision 
frequency  of  6  x  l(r  sec"*-  for  ground  state  CO  at  the  partial  pressures 
of  He,  CO,  and  sir  used  and  a  temperature  of  90°K.H  On  the  other 
hand,  assuming  that  all  the  electrical  energy  into  the  discharge  is 
dissipated  into  heat  which  is  distributed  uniformly  across  the  laser 
tube,  we  estimate  a  temperature  of  94°K  at  the  center  of  the  tube. 


Once  the  lineshape  function  was  determined,  a  second  experiment 
was  performed  to  get  more  accurate  values  for  the  parameters  |Ry_jJ  Ny 
and  |R^  Ny  ..  Oscillation  ranges  for  all  the  lasing  transitions 
withinVeach  band  were  measured  at  a  large  number  of  cavity  losses. 
Relative  line  center  gains  at  each  L  +  t  were  calculated  using  Eq.  (4) 
with  the  previously  determined  values  of  Av  and  T.  Relative  gains 
obtained  at  different  cavity  losses  were  al?  normalized  to  one 
arbitrary  scale  so  that  averages  could  be  taken.  Two  runs,  each 
covering  four  consecutive  bands,  were  made  under  slightly  different 
gain  medium  conditions.  The  respective  averaged  results  are  summarized 
in  two  columns  in  Table  1.  The  numbers  given  for  each  run  are  on  the 
same  arbitrary  scale,  but  the  results  for  the  two  different  runs  should 
not  be  compared  directly. 


The  linecenter  gains  given  in  Table  1  were  least  square  fitted  by 
Eq.  (1)  without  the  lineshape  factor  to  yield  |Ry_]J2  Ny  and  |Ry  ,  |2  Nv_i. 
The  least  square  error  was  calculated  according  to  Eq.  (5)  with  the 

temperature  T  varied  l°K1aS  a  time.  The  minimum  in  occurred  at  the 

same  temperature  for  all  bands  to  within  a  few  degrees.  The  T  giving 
the  smallest  y2  regardless  of  band  was  chosen  as  the  true  temperature 
for  each  run.  ^rfiey  were  T  =  91°K  and  93°  K  for  the  first  and  second 
runs  respectively.  The  experimentally  measured  gain  distribution  and 
the  best  theoretical  fit  for  the  v  =  7  -*  v  =  6  band  are  shown  as  the 
points  and  the  solid  curve  in  Fig.  4.  The  largest  deviation  of  an 
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experimental  point  from  theory  is  about  1%.  Finally,  ratios  of  the 
best  |rv+1 | 2  and  |r^_  y-  at  the  true  temperature  were  taken,  and 

we  obtained  the  ratios  of  the  matrix  elements  shown  as  open  circles 
in  Fig.  5.  Each  circle  in  this  figure  actually  represents  a  grand 
average  inasmuch  as  it  is  derived  from  a  statistical  fit  of  the 
relative  rotational  gain  distributions  of  two  consecutive  bands  which 
are  averages  of  data  obtained  at  different  cavity  losses. 

V.  ANALYSIS 

The  matrix  elements  in  Eq.  (1)  are  defined  by 

C1  */v<r>  Vi dr  «> 

where  Y  and  Y  -  are  the  vibrational  wave  functions,  and  the  dipole 
moment  function  H(r)  is  usually  expressed  in  a  Taylor  series  expansion 
about  the  equilibrium  nuclear  separation  r  as 


M(r)  Mt  (r  -  r^1  .  (7) 

i*0 

For  the  CO  molecule,  the  coefficients  have  been  determined  through 
i  =  3  from  absorption  intensity  data  onthe  fundamental  and  the^. 
first  two  overtone  bands.  >12  However,  the  matrix  elements  III1 12 
calculated  by  Young  and  Eachus  using  their  derived  dipole  moment 
coefficients  are  noticeably  smaller  than  those  implied  by  our  results. 
The  discrepancy  is  perhaps  to  be  expected  in  view  of  the  large 
uncertainties  in  the  M  .  Moreover,  probably  the  quartic  term  in 
2q.  (7)  should  be  included  in  the  evaluation  of  Ir^*  p  for  the  high 
v  levels  needed  in  the  present  work.  7 

An  alternative  way  of  expressing  the  dipole  moment  function  is  to 
expand  it  in  terms  of  all  the  matrix  elements  connecting  the  ground 
vibrational  level  as  follows:  ’ 


M(r)  Ro  vv 


(8) 


i=0 


The  matrix  elements  for  the  fundamental  bands  can  then  be  written  as 


(9) 
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Because  of  thi;  rapid  decrease  of  RQ  with  i,  only  the  first  few  tepns 
are  needed  in  the  actual  calculation.  After  Trischka  and  Salwen,  we 
shall  refer  to  this  approach  as  the  wave  function  approximation.  The 


iti.  1 

integrals  (i) 


evaluated  by  Herman  and  Rubin 
They  obtained  the  formula 


^  \dr  appearing  in  Eq.  (9)  have  been 

.  n  1  I  14  £■ .  1 _£  ..1 If _ * _ J  1  1  _  J 


for  the  case  of  the  Morse  oscillator. 


V(S-2i-l) 
i!  (v-1)!  v! 


2i  (S-2vf2)y_1  (S-2v)y 


(S-2i)^  (S-2v-i)i+1(S-2v-i)i_1 


where 


v-1  v 

ZE 

<k=0  m=0 


(-v+1).  (-v)  (S-2v-i) . 

in  _ yi  m  / «  *  *  v 

m!  (S-2V+2).  (S-2v)  (2v"i-'t'”m)i 

■C  m 


V. 


(10) 


(a). 


a(afl)  (a+n-1) ;  n  ^  1 


i; 


n  =  0 


and 


S  =  4  D  /u>  , 
e  e 


where  D  is  the  dissociation  energy  referred  to  the  minimum  of  the 
potentill  curve  and  U)  is  the  vibrational  frequency  for  infinitesimal 
amplitude. 

We  chose  to  fit  our  results  with  the  wave  function  approximation 
instead  of  the  Taylor  series  expansion  for  two  reasons.  First,  the 
matrix  elements  R*  (i  -  1,  2,  3)  are  known  with  higher  accuracy  than 
the  coefficients  M.(i  ■  1,  2,  3).  Hence,  the  chances  of  Improving  the 
expansion  in  the  former  case  seemed  better.  Also,  the  task  of  evaluating 
the  coefficients  Cv+1(i)  is  simpler  than  doing  the  necessary  integrals 
for  the  Taylor  serYes  expansion. 

Since  only  the  first  three  Rq  in  Eq.  (9)  are  known,  an  approxima¬ 
tion  for  R^"*  of  the  following  form  was  used: 

C1  -Z  Ro  C1  (i)  +  Ro4*  C1  <4)-  (11) 

i=l 

The  notation  Rq  is  used  here  to  signify  that  the  effect  of  the 
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truncation  is  absorbed  by  the  last  term. 


Toth  et  ril  have  compiled  in  their  naner  , 
measured  Inteu.ltle.  of  tl  CO  i-0?  °he° 

mean  value  for  each  band,  we  calculated  the  ratios  IX  2/R  1  -  nnfii  ^ 

VV  -11  three  m.«U  el^«2  h'eve  been 

£“  be  the  same  in  Ca8hion' s  convention13  by  Young  and  Eachus  12 
These  ratios  were  used  in  the  data  analysis. 


u_„  nf16  coefficients  cJfl(i)  were  evaluated  with  S  =  4D  /ou  =16 

with  sufficient"*  ^  necessary  Cv  (i>  could  be  calculated 
25  ?  «CCUraCy  ^6n  d°Uble  Praci8l«n  arithmetic  retaining 

n  Snificant  figures  was  employed.  The  ratios  |r^|2/|rv  |2  w  r 

~rtha  ya  ia8t  8qre  fJtted  t0  the  exPerimental  da?awi t^the 
i?  t  experimental  points  have  equal  uncertainties 

The  smallest  mean  square  error  occurred  for  rJ**/R  1  =  0  00014  in  the 

fZ  p°nVn^°n  °f  Cashion2  ratios  of  ma?rix  Siemens  catenated 

X2Pb(  in  the 

are  shin  e.  dot.  **£>  £  S^T^'i 

is  such  that  the  contribution  of  higher  R$  to  R^  isVgreater  for 

experlmenta^point^with  Mghi^  °^2  e““  f  *  ‘ha 

figure?  the  b6“  C‘1C“Uted  "“1“  -d  expert  t? 


167^  by 


For  v  -  10  each  term  in  Eq.  (11)  is  roughly  one  half  of  its 

I"  “*”ltude  has  the  opposite  sign.  If  thu  trend 

determination  of*,*  tor  the  CO  ^S?" 

.  the  opening  statement  of  the  last  paragraph,  one  sees 

immediately  that  a  10%  change  in  r^*/r  1  „dii  uZ  P  V f  L 

in  r  2/d  1  ro,  .  1?,K01  /K0  will  be  caused  by  a  2%  chance 

t  °<  5%  Change  in  “oW-  On  the  other  hand,  a  1%  change 

in  a^°f.the  efPeria*?tai1  ratios  in  the  same  direction  would  result 

/R°  •  “  th“«  aPPears  that  the  principal 

w  r  s\  may  come  from  experimental  errors  in  the  present 

nnint  5  ahows  that  the  largest  deviation  of  an  experimental 

present T  ?  u  ^  **  ab°Ut  2%*  The  laaer  tube  nsed  in  the 

srSiiiPw  ii  ::v5\at  its  *udut~  ^  *  • 

«li!h?  \  i  77  r  *  1  therefore  expected  that  there  would  be 

«  well  asainVth  ^  *“  ***  vibrational  level  population  densities 

as  well  as  in  the  temperature.  A  laser  tube  with  a  more  uniform 

active  medium  should  yield  improved  results. 

1 

Once  RQ  /R  i8  known,  the  matrix  elements  |rV  I2  can  be 

norasli  a  :ithM  aCCUracy  from  E(l-  dl).  The  results  when 
rmalized  to  (Rq1!  are  displayed  in  Table  2.  The  variation  of 
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|r'  |2  with  v  la  seen  to  be  almost  exactly  linear,  indicating  that  the 

effect)  of  mechanical  and  electrical  anharmonlcltles  almost  exactly 
cancel  each  other.  Also  shown  in  the  table  are  the  normalized  Einstein 
A  coefficients.  Due  to  the  frequency  factor,  they  increase  much  less 
rapidly  with  v. 


VI.  CONCLUSION 


In  summary,  we  have  presented  a  new  method  for  measuring  the 
relative  transition  probabilities  ol'  two  consecutive  molecular  bands. 
The  principal  requirement  for  good  accuracy  of  the  technique  is  that 
all  three  vibrational  levels  involved  be  nearly  equally  populated. 

An  application  of  this  method  to  the  fundamental  bands  in  the  CO  ground 
electronic  state  has  yielded  for  the  first  time  the  ratios  Ir7*  |2/ 

|RV  1 |2  for  seven  consecutive  bands.  The  experimental  ratios 
substantiate  earlier  calculations  of  |r^  , | 2  for  the  same  transitions. 
They  imply  an  almost  exactly  linear  dependence  of  ]R^  -ll  on  v* 


It  is  of  Interest  to  compare  the  remits  from  the  present  work  with 
the  theoretical  calculations  of  Cashion.  By  equating  Riv  from  the 
harmonic  oscillator  model  to  R^v  from  the  Morse  oscillator  model,  he 
obtained  an  expression  for  Rqv/RqV  .  He  showed  that  all  the  Rqv 
should  have  the  same  sign  for  sufficiently  low  v.  This  has  been 
demonstrated  to  be  the  case  for  CO.  as  we  mentioned  earlier,  for  the 
first  three  v  by  Young  and  Eachus.  2  The  present  experiment  extends 
Cashion' s  contention  to  v  *  4.  The  magnitude  of  Rq^  obtained  by  us  is 
only  about  2/3  of  the  theoretical  value  calculated  by  Cashion.  However 
as  we  have  indicated  in  the  last  section,  our  value  for  Rq**  could  be 
in  error  by  a  fairly  large  fraction. 


The  method  of  measuring  relative  transition  probabilities  described 
inln  this  paper  can  be  applied  most  immediately  to  the  other  existing 
molecular  laser  bands.  In  cases  where  the  single  pass  optical  gain  is 
large,  it  may  be  more  convenient  to  use  the  direct  probe  method  for 
measuring  gain  rather  than  the  laser  oscillation  range  method  which 
Involves  a  fair  amount  of  work  in  data  reduction. 


It  is  worthwhile  to  point  out  that  an  elaborate  experimental 
procedure  was  adopted  In  the  present  work  only  because  we  desired  an 
accuracy  close  to  one  percent  for  the  ratios  of  the  matrix  elements. 
For  many  applications  considerably  lower  accuracies  would  probably 
suffice. 


The  approach  presented  in  Section  II  can  also  be  used  to  deduce 
relative  matrix  elements  for  parallel  bands  which  share  the  same  upper 
or  lower  level.  We  emphasize  that  the  method  is  equally  applicable  to 
molecular  bands  which  are  mostly  absorbing.  Thus,  relative  transition 
probabilities  in  molecular  systems  with  optical  loss  can  also  be 
measured.  This  is  an  important  point  for  the  widespread  use  of  the 
method,  since  only  a  limited  number  of  molecular  bands  have  been 
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observed  to  exhibit  optical  gain.  Experiments  on  non- lasing  bands 
should  be  practicable  with  the  use  of  a  tunable  IR  laser.  Investigations 

of  this  type  should  prove  to  be  particularly  interesting  for  polyatomic 
molecules. 
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FIGURE  1 

Illustration  of  single  mode  laser  oscillation  condition 
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FIGURE  2 

Illustration  of  laser  oscillation  ranges  for  a  distribution 
of  vibration-rotation  transitions. 
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FIGURE  5 
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Ratio  of  matrix  elements  for  adjacent  CO  fundamental 
bands  as  a  function  of  v.  Open  circles  are  experimental 
values;  dots  are  best  theoretical  fit. 
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TABLE  1 


Experimental  llnecenter  gain  distributions.  Results  in 
the  two  columns  are  from  two  separate  runs.  Each  set 
is  on  an  arbitrary  scale. 


V+l-V 


12-11 


11-10 


10-9 


9-8 


P(J) 

exp 

_ °P(J) 

V+l-V 

PCI) 

exp 

°P(J) 

P(5) 

0.890 

9-8 

P(5) 

1.224 

P(6) 

1.325 

P(6) 

1.834 

P(7) 

1.483 

P(7) 

2.066 

P(8) 

1.486 

P(8) 

2.094 

P(9) 

1.235 

P(9) 

1.803 

P(10) 

0.950 

P(10) 

1.425 

P(ll) 

0.669 

P(ll) 

1.050 

P(5) 

0.995 

8-7 

P(5) 

1.227 

P(6) 

1.499 

P(6) 

1.897 

P(7) 

1.740 

P  ( 7) 

2.230 

P(8) 

1.669 

P(8) 

2.215 

P(9) 

1.364 

P(9) 

1.902 

P(10) 

0.942 

P(10) 

1.481 

P(ll) 

0.712 

P(6) 

1.612 

7-6 

P(5) 

0.908 

P(7) 

1.904 

P(6) 

1.660 

P(8) 

1.833 

P(7) 

2.060 

P(9) 

1.517 

P(8) 

2.092 

P(10) 

1.060 

P(9) 

1.869 

P(ll) 

0.705 

P(10) 

1.465 

P(ll) 

1.047 

P(5) 

1.053 

6—5 

P(7) 

1.318 

P(6) 

1.725 

P(8) 

1.499 

P(7) 

1.995 

P(9) 

1.509 

P(8) 

1.991 

P(10) 

1.266 

P(9) 

1.637 

P(ll) 

0.945 

P(10) 

1.165 

P(ll) 

0.792 
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TABLE  2 


Matrix  elements  and  corresponding  Einstein  A  coefficients 
calculated  from  the  wave  function  approximation  using 
R^*  obtained  from  experimental  results. 


V 

K.ii2/i*o1i2 

AV  /A 

Vro 

1 

1.00 

1.00 

2 

2.00 

1.93 

3 

2.99 

2.78 

4 

3.98 

3.56 

5 

4.97 

4.28 

6 

5.96 

4.92 

7 

6.95 

5.53 

8 

7.95 

6.09 

9 

8.95 

6.59 

10 

9.97 

7.04 

11 

11.00 

7.46 

12 

12.00 

7.81 
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COLD  CATHODE  ELECTRON  BEAM  CONTROLLED  I.ASER  AMPLIFIER*^ 
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J.  Jenson  and  A.C.  Kolb 
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San  Diego,  California  92123 


ABSTRACT 


A  large  area  cold  cathode  (10  cm  x  100  cm)  electron  gun  has 
been  constructed  and  used  to  control  a  CO  laser  amplifier.  Energy 
deposition  of  226  J/ liter  in  a  time  of  1  /nsec  has  led  to  a  gain  of 
4. 97J cm  in  an  atmospheric  pressure  COj  laser  mixture. 


t  To  be  published  in  Applied  Physics  Letters. 
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.  k  en6nL  experi“ents  l,y  Fenstermacher,  Nutter,  Leland  and  Boyer' 
and  by  Daugherty  Pugh,  and  Douglas-Hamilton^2)  have  shown  that  an 

electro"  beam  atabilize  larSe  volume  C02  laser  discharges.  The 

electron  guns  for  these  experiments  used  thermionic  emission  from  hot 

a/ioder  or  tetrode^2)  structure.  Another  possiMe 
source  of  electrons  for  laser  control  is  the  cold  cathode.  In  contrast 
to  the  thermionic  cathode  which  is  capable  of  supplying  moderate 

cathod  den3ities  (<  -1  A/ cm2)  for  arbitrary  oulse  length  the  rold 

cathode  operates  satisfactorily  at  high  current  densities  (>  1  A/ cm2) 
or  me  periods  less  than  a  few  microseconds.  A  cold  cathode 
electron  gun  is  very  simple  in  design,  can  be  economical  in  certain 
PP  ications,  and  is  ideally  suited  for  pumping  laser  amplifiers  for 

“  °f  °Peratl0n-  ”-us  it  is  important  to  show  rtst 

??2  laser  mad*a  can  be  satisfactorily  pumped  under  control  of  these 
high  current  density,  short  duration  e-beams.  We  have  constructed  * 
large  area  (10  cm  x  100  cm)  planar  cold  cathode  electron  gun  and  Lve 
conducted  preliminary  experiments  on  the  control  of  a  discharge  in  a 

Uslng  3  3:^:1  mixture  of  He:N2:C02  at 

a oil ^cmTSST  W%haVe  mea8Ured  Sn  unsaturated  gain  coefficient 
a  4.9/o  cm  at  10.6  n,  for  an  energy  input  of  226  J/ liter  into  the 

cathodes\lO  iff  has  demonstrated  that  large  area  cold 

cathodes  (10  ^  x  10  cm)  can  supply  uniform  electron  current  densities 
of  1  -  10  A/u*2  at  200  keV  for  time  durations  of  up  to  5  Msec  Sese 

risinfvolt6/11  PilnCiple  V6ry  simPle*  0ne  “«ely  applies  a  ^apidly- 
rising  voltage  pulse  to  a  conductor  which  has  a  sharp  edge  facing  the 

anode.  The  most  likely  mechanism  for  diode  operation  begins  ^h 

wf  d  emission  fr°®  the  regions  of  high  field  stress  (>  10?  V/cm) 

centeyed  on  sharP  ed8es  or  whiskers.  The  rapid  heating  of 
these  microscopic  projections  leads  to  their  vaporization  and  the 
creation  of  a  conducting  plasma  which  serves  as  ?he  source  £  electrons 
Sen  Pjasma'  however,  causes  the  impedance  of  the  diode  to  collapse 

Before  theaSn8d  “d  Sh°rt  circuits  the  anode-cathode  space, 

limited  iL^I)  r  aollaPses  the  dioda  operates  in  a  space  charge 
i  f  m  f  and  tbe  r?}ftion  approximately  follows  Child's 

spacing  i/ cTZ/vT’.l  ""  *«■  d  is  the  cathode-anL 

the  current™  ^peree  ^  the  dl°d<i  ‘n  V°Us'  OTd  1 

hladp0^.-10^  V00  Cm  C°ld  cathode  was  constructed  from  strip  razor 

£? Vcm  w'hl t0  6nhanCe  Plasma  f°™ati°"  at  the  edge. 

Sd  S  V  8  bJades  were  oriented  transverse  to  the  axis  of  the  gun 

and  °n  fcfentfra-  7116  blades  were  mounded  on  the  corntrs 

plane  5  5  ^fr^theVn ' ^  h°lder  which  was  mounted  in  the  horizontal 
plane  5.5  cm  from  the  foil  support  structure  (see  Fig.  1)  on  speciallv 

designed  Teflon  insulators.  Hie  diode  was  connected  [o  a  six-stage 

Marx  generator,  which  erected  to  209  kV  in  a  time  of  100  nn  and  i,8d 

output  capacity  of  0.05  „F.  The  voUage  ^  in  1°  “ 

cathode  inductance  reduced  the  catSlde  voltage  by  44  kV.  Thia  yielded. 


(1) 
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after  0.3  tisec,  a  current  of  4  A/ cm  which  agrees  with  the  calculated 
space-charge  limit.  If  the  transmission  factor  of  the  foil  (.0015"  A-C) 
and  support  structure  are  taken  into  account ,  one  would  observe  a 
peak  current  density  of  ~  3  A/ cm2  in  the  laser  gas.  The  experimental 
current  density  is  estimated  to  be  0.3  to  0.8  A/ cm2  from  the  gas 
conductivity'*''  and  from  Faraday  cup  measurements.  The  deviation  from 
simple  geometric  transmission  through  the  foil  support  structure  is 
most  likely  caused  by  large  transverse  moments  imparted  to  the 
electrons  upon  emission  from  the  razor  blade  edges.  This  is  supported 
by  the  results  of  J.C.  Martin(),  which  indicate  that  the  half  angle 
of  electron  emission  from  a  razor  blade  cathode  is  about  45° . 


Electron  beam  uniformity  measurements  are  required  to  complete 
the  studies  of  these  electron  guns.  Moreover ,  higher  diode  voltages 
should  be  used  to  enhance  the  cathode  plasma  formation  process  and 
ensure  good  uniformity  of  emission.  Measurements  on  the  impedance 
collapse  of  the  diode  indicated  that  the  plasma  closed  toward  the 
anode  at  a  rate  of  about  4  cm//isec.  This  is  greater  than  the  rate 
measured  by  Loda  and  DeHart  which  is  ~  1.5  cm/ /nsec. 

The  laser  plasma  was  generated  at  atmospheric  pressure  between 
a  10  cm  x  100  cm  kt  anode  and  stainless  steel  wire  cathodes  (see  Fig.  1) 
The  electrode  spacing  for  most  experiments  was  5.0  cm.  The  gas  was 
contained  in  a  plexiglass  box  which  had  been  lined  with  absorbing 
material  on  the  inside  surface  to  prevent  parasitic  oscillation.  In 
addition,  the  NaCl  windows  were  canted  in  opposite  directions  and  the 
kt  anode  was  sandblasted  to  prevent  specular  reflection.  A  sealed-off 
probe- laser  utilizing  an  invar  cavity  was  operated  at  10.6  n  [P(18) 
to  P(22)]  to  probe  the  amplifying  medium.  The  probe  beam  was  directed 
through  the  amplifier  box,  into  a  screen  room,  and  on  to  an  Au:Ge 
detector  that  was  cooled  to  77°K.  The  detector,  bias  circuit,  and 
oscilloscope  system  had  a  risetime  faster  than  1  n sec,  and  the  system 
was  measured  to  be  linear  over  the  dynamic  range  of  the  gain  measure¬ 
ments.  The  laser  frequency  was  brought  to  line  center  with  a  piezo¬ 
electric  translator;  no  beam  wander  due  to  the  e-beam  produced  plasma 
was  evident. 


Figure  2  shows  the  gain  coefficient  plotted  as  a  function  of 
electrical  energy  deposited  into  the  gas.  In  the  insert  of  Figure  2 
is  a  typical  oscillogram  showing  the  time  history  of  the  sustainer 
current  pulse  and  development  of  the  optical  gains  Notice  that  the 
gain  peaks  in  the  afterglow  of  the  sustainer  discharge.  Variations  in 
the  energy  deposition  into  the  laser  plasma  were  achieved  by  changing 
the  E/P  ratio  for  a  fixed  pumping  time  of  typically  ~  1  /isec.  The 
peak  gain  was  measured  for  an  E/P  ratio  of  6.6  kV/cm.  Extrapolating 
from  Fenstermacher '  s  data''-',  a  0.5  A/on^  electron  beam  will  control 
about  39  A/ cm2  Qf  sustainer  current  (E/P  5.0  kV/atm-cm).  This  is 
in  substantial  agreement  with  our  results.  Under  these  conditions, 
sustainer  energy  is  delivered  to  the  gas  at  the  rate  of  190  J/liter//isec 
Limitations  of  charge  storage  in  the  Marx  bank  driving  the  diode 
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precluded  investigations  where  the  conductivity  of  the  plasma  was 
maintained  for  times  greater  than  1  /nsec.  The  optical  energy- storage 
efficiency  as  reflected  by  the  gain  is  lower  than  that  observed  by 
Fenstermacher  et  al.u'  This  may  be  due  to  the  different  rates  of 
energy  deposition  in  the  two  devices.  The  scatter  in  the  data  is  most 
likely  due  to  a  reduction  in  the  e-beam  range  in  the  gas  as  the  cathode 
voltage  dropped  and  to  random  nonuniformities  in  the  e-beam  as  the 
diode  impedance  collapsed.  Since  the  probe  laser  sampled  only  a  very 
small  volume,  any  deviation  in  the  pumping  rate  over  this  volume  shows 
up  in  the  total  small  signal  gain. 

The  results  of  these  experiments  indicate  that  the  cold-cathode 
electron  gun  should  offer  a  relatively  simple  means  of  controlling 
large-volume  gas  lasers.  This  method  is  particularly  useful  in 
situations  where  rapid,  high  current  pumping  is  desirable  «2  -  3  /xsec), 
as  is  the  case  in  COo  amplifiers  designed  for  nanosecond  pulse 
amplification.  At  the  measured  gain  of  4.9%  cm-1,  5  to  6  joules  of 
optical,  10.6  fi  energy  is  stored  in  the  medium  which  can  be  extracted 
in  a  nanosecond  pulse. 
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SUSTAINER 


FIGURE  1.  Schematic  diagram  of  experimental  apparatus 
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APPENDIX  C 


DAMAGE  MEASUREMENTS  WITH  SUBNANOSECOND  PULSES* 

John  M.  McMahon 

Naval  Research  Laboratory 
Washington,  D.C.  20390 


1.  Introduction 

There  is  currently  a  great  deal  of  interest  in  generating 
intense  subnanosecond  pulses  for  laser  CfR  applications  as  well  as 
for  other  purposes.  Success  in  construction  of  relatively  efficient 
solid  state  lasers  which  generate  terawatt  (or  higher)  level  pulses 
and  which  operate  reliably  for  many  shots  will  require  either  that  the 
pulse  shape  and  intensity  be  tailored  in  time  and  spatial  extent  such 
that  non-linear  effects  do  not  predominate  or  that  geometries  be 
chosen  in  which  these  effects  can  reasonably  be  expected  not  to  occur. 

Disc  amplifiers,  for  example,  appear  to  be  relatively  immune  to 
self  focusing  (if  only  in  that  the  surfaces  will  damage  first).  They 
however  are  somewhat  unattractive  from  an  efficiency  and  convenience 
standpoint  compared  to  a  solid  rod  system. 

We  have  investigated  the  self  focusing  problem  for  three  laser 
glasses  in  common  use  (Owens-Illinois  ED-2J  Soveril  MG  915  and  Schott 
LG-56) .  Measurements  of  N2  are  reported  on  unpumped  rods  for  pulse 
durations  between  20  psec  and  250  psec.  Measurement  techniques  for 
subnanosecond  pulses  are  not  presently  as  quantitative  as  might  be 
desired.  Pulse  duration  measurements  by  non-linear  optical  techniques 
may  lead  to  erroneous  results  unless  great  care  is  taken.  Energy 
measurements  are  complicated  by  plasma  formation  in  the  calorimeter  at 
very  low  levels  compared  to  q  switched  pulses  (^  50  mJ/cm  on  carbon) 
which  reduces  the  signal  to  noise  ratio  to  (typically)  less  than  20:1 
and  there  are  difficulties  in  measuring  near  and  far  field  intensity 
distributions  with  a  high  precision  and  wide  dynamic  range  at  1  jan. 

2.  Measurement  of  No  in  Laser  Glasses 

In  the  1970  ASTM  symposium,  J.  Davit^  of  CGE  and  M.  Duguay  and 
J.  Hansen(2)  of  BTL  reported  on  measurements  on  N2  for  several  glasses 
by  various  methods.  Duguay  and  Hansen  reported  some  measurements  at 
~  8  psec  on  BK-7  and  LASF-7  glasses  and  Davit  on  tests  done  at  2  nsec 
on  LG-56  and  MG-915  glasses.  While  BK-7  might  be  expected  to  behave 
in  a  similar  fashion  to  the  laser  glasses  Davit  tested,  the  value  of 
N2  was  about  an  order  of  magnitude  larger  than  the  BTL  value.  This 
suggests  the  possibility  that  for  very  short  pulses  N2  might  be 

*This  work  was  supported  by  the  Advanced  Research  Projects  Agency 
under  ARPA  Order  2062. 

Paper  presented  at  the  4th  ASTM-NBS  Symposium  on  Damage  in  Laser 
Materials,  Boulder,  Colorado,  June  14-15,  1972. 
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smaller  due  to  mechanisms  Important  on  the  nanosecond  time  scale  not 
being  able  to  respond  to  short  pulses.  Another  possibility  was  that 
since  neither  author  used  a  laser  which  produced  a  time  band  width 
limited  pulse,  both  measurements  could  be  consistent  with  a  model  in 
which  some  of  the  pulse  energy  is  in  very  short  temporal  spikes. 

For  the  measurements  reported  here  we  used  a  mode  locked  Nd:YAG 
oscillator-amplifier  which  has  been  described  elsewhere. (3)  This 
system  produces  a  TEMQO  spatial  mode  output  single  pulse  with  an 
energy  adjustable  from  1-200  mJ  in  a  pulse  whose  width  can  be  adjusted 
from  20  psec  to  1.3  nanoseconds  by  the  use  of  transmission  etalons  in 
the  cavity.  Our  measurements  and  those  of  Carmen  and  co-workers  at 
LLL  indicate  that  the  pulse  is  essentially  time-bandwidth  limited™) 
with  a  smooth  Gaussian  spectral  and  temporal  output,.  Figure  1  shows 
the  geometrical  arrangement  of  this  device. 


For  the  trapping  measurements  the  bea®  was  propagated  ~  3  meters 
to  get  rid  of  near  field  Fresnel  fringes  and  the  pattern  checked  and 
found  to  be  a  Gaussian  down  to  below  the  e-3  points.  The  beam  was 
then  run  through  a  sample  rod  and  the  trapping  length  noted  from  the 
side  during  10-20  shots  for  each  sample.  In  all  cases  the  start  of 
the  filament  was  characterized  by  the  emission  of  blue  light  (presum¬ 
ably  from  plasma  breakdown  in  the  track).  This  was  somewhat  enhanced 
in  the  ED-2,  presumably  due  to  the  ultraviolet  pumping  the  cerium  which 
then  fluoresced  in  the  region  of  4000-5000  A. 


Table  1  lists  the  geometries  and  absorption  coefficients  of  the 
samples  tested.  The  Soveril  glass  tested  was  a  rod  made  by  the  old 
process  and  had  a  significantly  higher  absorption  coefficient  than 
MG  915  made  by  the  continuous  flow  process.  The  LG  56  rod  was 
supplied  by  CGE  with  the  NRL  VD  640  and  had  a  low  absorption 
coefficient  for  LG-56.  We  initially  questioned  whether  this  was  LG  56 
or  LG  630  and  were  assured  by  CGE  that  it  was  LG  56. 

In  the  actual  sequence  of  experiments  the  laser  was  tuned  to 
fUn  mlon  ®^ven  Pulsewidth  and  then  the  samples  were  each  irradiated 
tor  10-20  shots  and  the  mean  trapping  length,  peak  energy  density, 
pulsewidth  and  spatial  distribution  determined.  The  laser  was  then 
tuned  to  the  next  pulsewidth  and  the  sequence  repeated.  Table  2  shows 
the  results  for  puls©  durations  of  20  psec  and  250  psec.  The  pulse 
peak  powers  and  beam  parameters  were  PQ  =  8  x  1010  W/cm2,  b  -  1.1  am 
and  PQ  “  1  x  10  W/cm*,  b  **  1.5  mm,  respectively. 

The  difference  in  beam  parameter  was  caused  by  the  beginning  of 
focusing  in  the  YAG  preamplifiers  which  was  not  severe  enough  to 
noticably  perturb  the  temporal  or  spectral  structure  of  the  pulse.  The 
cone  angle  of  the  self-phase  modulated  light  from  the  samples  was 
““  0.1  radian  at  250  psec  and  0.25  radian  at  20  psec.  Typically 
3-10  filaments  would  be  formed  which  would  traverse  the  length  of  the 
rod. 
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Table  3  gives  trapping  data  on  ED-2  for  a  somewhat  more  complete 
set  of  experiments.  Data  sets  were  run  for  20  psec,  125  psec  and 
250  psec  "clean"  pulses  and  also  for  a  pulse  125  psec  long  with  the 
etalon  misaligned  to  give  two  spectral  components  of  equal  intensity. 
These  interfered  giving  a  peak  intensity  ~  four  times  the  intensity 
for  a  non- structured  pulse.  This  caused  a  reduction  in  trapping 
length  of  a  factor  of  two. 

Table  4  gives  the  calculated  value  of  N2  for  the  glasses  tested 
with  an  estimated  absolute  accuracy  of  1407..  The  value  of  BK-7 
should  be  taken  as  a  upper  limit  since  -68ampie  Zfoc 

3 .  Discussion  of  the  Results 

The  data  all  showed  an  excellent  agreement  with  Kelley's  formula 
for  the  parametric  variation  of  self  focusing  length^ 


Additionally,  over  the  range  of  temporal  values  tested  N2  for  a 
particular  glass  was  found  to  be  sensibly  constant. 

It  was  additionally  shown  that  N2  will  follow  temporal  structure 
on  the  pulse.  If  we  had  assumed  our  structured  pulse  to  be  smooth,  an 
N2  value  four  times  the  observed  value  would  have  been  obtained. 

Under  a  given  set  of  conditions  there  did  appear  to  be  small  but 
experimentally  significant  differences  in  self  focusing  length.  That 
there  is  a  difference  between  ED-2  and  LG-56  is  not  too  surprising 
since  the  two  are  dissimilar  glasses,  but  the  difference  between  the 
LG-56  and  MG-915  is  surprising,  since  they  are  essentially  the  same 
glass.  These  may  be  related  to  the  difference  in  absorption 
coefficient  of  the  two  glasses. 

4.  Systems  Considerations 

Based  upon  the  values  of  N-  reported  here  there  appear  to  be  no 
fundamental  difficulties  in  amplifying  pulses  to  high  energies  if  the 
trapping  length  is  long  compared  to  the  amplifier  length  and  the 
pulse  is  temporally  and  spatially  smooth.  Unfortunately,  these 
conditions  are  rarely  met  (for  more  than  a  few  shots). 

Dust  particles  on  surfaces  can  provide  nuclei  for  wavefront 
disturbances  which  can  then  cause  trapping  in  the  next  stage.  Shot  by 
shot  this  damage  can  then  propagate  through  the  system.  This  problem 
can  be  minimized  if  the  amplifiers  are  spaced  in  the  far  field  of  any 
filaments  generated  in  the  previous  amplifier. 


C-3 


TABLE  1 


SAMPLES  TESTED 


GLASS  TYPE 


ED-2 

MG-915  (cast) 

LG-56 

BK-7 


4Acm)  V(cm^) 


30 

30 

50 

6 


3  x  lO-3 
7.5  x  10 
3  x  10", 
3  x  10 


TABLE  2 


RELATIVE  DATA  ON  TRAPPING  LENGTH 


TYPE 

Z(250  psec) 

Z(20  psec) 

ED-2 

24  ±  1 

6.6  ±  1 

MG-915 

24  ±  1.5 

6.1  ±  1 

LG- 56 

21  ±  1.5 

5.6  ±  1 

BK-7 

6  +  * 

*  The  6  cm  sample 
within  1  cm  of  the 

only  trapped  occasionally  (=~ 
surface. 

TABLE  3 

TRAPPING  DATA  ED-2 

307.  of  shots) 

t(psec) 

PpMcm  2) 

b(mm) 

Z(cm) 

20 

8  x  1010 

1.1 

6.6  ±  . 

125 

1.2  x  1010 

1.2 

18  ±  2 

250 

1  x  1010 

1.5 

24  ±  1 

125* 

4.8  x  1010 

1.2 

9.2  ±  1 

*  Etalon  misaligned  to  cause  temporal  beats  with  a  four- fold  enhance¬ 
ment  of  peak  power. 
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TABLE  4 


BK-7 

ED-2 

MG-915 

LG-56 


4 


VALUES  (eau) 


6  x  10 
(7  ±  3)  x  10 
(7  ±  3)  x  10 
(8  ±  3)  x  10 


-14 

-14 

-14 
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It  Is  also  generally  very  difficult  to  avoid  Fresnel  fringes  In 
the  beam  pattern.  Based  on  the  experiments  of  Gulliani,(6)  a  Fresnel 
fringe  Itself  may  be  expected  to  have  a  high  threshold  for  trapping 
(since  It  would  be  a  one  dimensional  process).  However ,  In  a 
circularly  symmetric  beam,  the  Fresnel  fringes  may  accentuate  damage 
by  moving  energy  Into  the  beam  center  where  it  may  tend  to  form  an 
Intensity  spike.  We  have. seen  such  a  spike  forming  at  times  and  we 
understand  that  Bettlnger'  '  at  CGE  has  seen  such  a  spike  form  and 
then  self-focus.  This  effect  seems  to  be  somewhat  diminished  If  any 
apertures  in  the  system  are  non- concentric. 

Operating  the  CGE  VD  640  amplifier  train  with  250  psec  pulses 
outputs  of  50-60  joules  have  been  achieved  with  little  or  no 
degradation  over  **  100  shots.  A  limited  number  of  shots  have  been  run 
at  levels  of  90-100  shots.  No  degradation  has  been  observed  to  date. 
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1.  INTRODUCTION 


Recent  studies  of  the  construction  of  large  laser  systems  have 

emis«lnn°ULthei?eed  cfreful  consideration  of  amplified  spontaneous 
emission  and  self-oscillation  within  the  laser  medium  itself.  These 

effects  can  seriously  reduce  the  efficiency  of  the  laser  system  by 
extracting  energy  from  the  gain  medium  in  undesirable  modes 
Proposals  for  iarge  Nd  doped  glass  laser  amplifiers  are  suggesting 
glass  disks  of  up  to  25  cm  in  major  diameter  and  up  to  several  cm 
thick.  Internal  "whisper"  modes  within  these  disks  and  amplification 
of  spontaneously  emitted  light  that  is  scattered  off  of  the  edges  or 

capabilit  yf  ^  di8kS  ^  8erlously  reduce  their  energy  storage 

? 6  have  aome  simple  measurements  on  the  internal  scattering 

TlliniiiC|In0J  f°U8h  8fOUnd'  Polished,  coated  surfaces  of  1%  Nd-Owens 
Illinois  ED-2  laser  glass.  These  measurements  have  been  made  using 
an  integrating  sphere  technique.  They  point  out  the  wide  variation 
in  scattering  efficiency  of  glass  surfaces  as  a  function  of  angle 
and  surface  preparation.  These  measurements  in  conjunction  with  recent 
computer  codes  that  can  trace  rays  in  arbitrary  geometries  with 
arbitrary  absorption  coefficients  (including  gain)  can  be  used  to 
estimate  the  extent  of  the  above  problems. 

2.  MEASUREMENT  TECHNIQUE 

„  FJ?U5e  1  ®hows  the  technique  that  was  used  to  measure  the  internal 
reflectivity  of  various  suriaces.  The  incoming  laser  beam  enters  the 
ma  integrating  sphere  (a  ping-pong  ball  painted  with  BaSO.  paint) 
and  scattere  off  of  both  glass  surfaces.  We  are  interested  only  in 
the  scattering  from  the  rear  surface.  The  front  surface  reflection  is 
allowed  to  escape  from  the  sphere  through  a  small  hole  that  is  suitably 
positioned  and  the  rear  surface  reflection  scatters  within  the  sphere. 
After  the  first  bounce  off  of  the  sphere,  the  interior  of  sphere  is 
uniformly  illuminated  and  subsequent  ray  trajectories  are  statistically 
identical.  This  is  independent  of  the  specular  or  diffuse  nature  of 

^  reflectivity-  The  germanium  photodiode  (Electro-Nuclear 

No.  653)  i8  pointed  so  that  only  scattered  light  is  received.  It 
detects  the  relative  amount  of  scatter  from  various  glass  samples. 

The  sphere  can  be  calibrated  by  using  a  BaS(>4  surface  to  provide 
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total  diffuse  scattering.  The  Incoming  laser  beam  Is  chopped  at 
200  Hz  with  a  tuning  fork  modulator  and  the  modulated  signal  from  the 
photodiode  is  detected  with  a  lock-in  amplifier.  A  1.15  micron  He-Ne 
laser  was  used  because  of  its  availability  and  because  most  scattering 
and  absorbing  processes  in  the  materials  that  we  have  examined  are 
essentially  independent  of  wavelength  in  the  1.06  to  1.15  micron 
region. 

The  linearity  of  the  system  was  verified  with  known  attenuators 
and  the  scattering  accuracy  was  verified  by  comparing  the  front 
surface  Fresnel  reflection  from  a  piece  of  ED-2  glass  to  a  BaSO^ 
covered  plate.  The  67.  reflectivity  agrees  with  the  Cary  14  measure¬ 
ment  of  5.57.. 

An  error  is  introduced  in  that  diffusely  scattered  radiation 
from  the  rear  surface  that  is  below  the  critical  angle  will  be 
trapped  within  the  glass.  The  extent  of  this  effect  can  only  be 
estimated  in  the  present  system  and  is  discussed  in  Section  3.  Care 
must  be  taken  to  ensure  that  only  laser  radiation  that  is  scattered 
from  the  back  surface  Js  measured.  A  slight  amount  of  scattering 
of  the  front  surface  reflected  beam  from  the  edge  of  the  exit  hole 
could  lead  to  erroneous  results.  Effects  to  reabsorption  of  scattered 
light  by  the  test  surface,  and  from  non-specular  scattering  from  the 
front  surface  have  not  been  considered.  These  effects  would  be  small. 

3.  EXPERIMENTAL  RESULTS 


The  data  presented  in  the  following  Table  (Table  1)  was  taken 
with  the  experimental  apparatus  that  was  discussed  in  Section  2. 

The  samples  were  0.63  cm  thick  slabs  of  Owens  Illinois  ED-2  glass 
that  were  smooth  ground  (22  n  AI5O3  abrasive)  one  side  and  smooth 
polished  (pell  polish)  on  the  otner.  Sample  No.  2  was  coated  with 
Owens  Illinois  black  glass  on  the  ground  surface.  Sample  No.  1  was 
not  coated  at  all,  and  Sample  No.  3  was  smooth  polished  on  both  sides 
and  then  coated  on  one  side  with  the  Owens  Illinois  black  glass.  On 
this  sample  there  was  no  diffuse  reflection.  The  incident  beam  was 
13°  from  the  normal.  The  reflectivities  are  corrected  for  the  front 
surface  Fresnal  losses. 

The  black  glass  coating  is  black  "solder"  glass  that  was  specially 
developed  for  the  purpose  of  cladding  ED-?  laser  glass.  It  has  no 
Owens  Illinois  designation  to  date. 
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TABLE  I 


SAMPLE  NO. 

REFL.  (|U  VOLTS) 
SIGNAL 

%  REFL.  COMPARED 
TO  BaS04 

%  REFL.  CORRECTED 
FOR  CRITICAL 

ANGLE  EFFECTS 

BaSO^ 

4200 

100% 

100% 

1.  Smooth  Grnd. 

No  coating 

100 

2.6% 

6.5% 

2.  Smooth  Grnd. 
0.1.  Coating 

45 

1.2% 

2.7% 

3.  Smooth  polish 
O.I.  Coating 

10 

0.27% 

0.27% 

Front  surface  refl 

.  220 

5.2% 

5.2% 

Smooth  grnd.  to  laser  390 
beam 

9.2% 

9.2% 

The  data  is  reproducible  to  within  ±  10%,  and  on  measurements 
that  can  be  checked  with  the  Cary  14,  the  agreement  is  good.  However, 
the  corrected  data  depends  on  the  surface  being  Lambertian.  Angular 
measurements  indicate  that  at  normal  incidence  this  may  be  reasonable, 
however,  this  data  is  certainly  less  accurate,  probably  ±  25%. 

The  amount  of  light  that  can  escape  from  a  Lambertian  surface 
inside  a  medium  of  higher  index  is  restricted  by  total  internal 
reflection  at  low  angles. 


If  the  light  is  conducted  out  of  the  measuring  region  and 
completely  reabsorbed  as  in  our  system,  the  detected  signal  must  be 
corrected  for  the  lost  light.  This  effect  is  non-negligible. 

For  a  Lambertian  surface  that  scatters  as  C0S%\ 


.  Bln2  6 

*in 


.41 ! 
'n=l. 


56 
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The  detected  signals  from  the  diffuse  absorber  are  corrected  by 
multiplying  column  3  by  2.4  and  are  placed  in  column  4  of  Table  I. 

Low  Angle  Scattering 

Additional  data  on  the  low  angle  scattering  efficiency  is 
presented  in  Table  II.  This  data  was  taken  by  comparing  the  reflected 
light  to  the  inching  laser  light  as  measured  directly  with  the 
integrating  sphere.  See  Figure  2.  Tht  laser  polarization  was  perpendi¬ 
cular  to  the  plane  of  incidence. 


TABLE  II 


Scatter  from  Following  Raw  Signal  Percent  Trans.  As 

Sample  m  Volts  Function  of  Internal 


Angle 

Signal 

Angle  From  Surface 
(Corrected  for  Fresnel 
Losses) 

Angle  Refl. 

Laser  Signal 

50 

1007. 

Attenuation  thru  2  surfaces 

44 

88% 

No.  3  Smooth  polish 

1.6° 

33 

1.6° 

747. 

spec. 

0.1.  black  glass 

3.2° 

25 

3.2° 

547. 

spec. 

6.4° 

17 

6.4° 

387. 

spec. 

No.  2  Smooth  grind 

9.5° 

12 

9.5° 

27% 

spec. 

3.2° 

21 

3.2° 

47% 

mostly 

0.1.  black  glass 

6.4° 

12 

6.4° 

277. 

spec. 

mostly 

Smooth  polish 

3.2° 

35 

3.2° 

797. 

spec. 

mostly 

ED-2  surface 

6.4° 

39 

6.4° 

85% 

spec. 

mostly 

spec. 

In  this  configuration,  light  that  is  scattered  at  large  angles 
will  not  be  detected  and  this  can  lead  to  some  errors.  The  data  in 
the  last  column  for  a  polished  ED-2  surface  should  be  1007.,  because 
the  beam  would  undergo  total  internal  reflection.  This  indicates 
that  some  light  has  not  been  accounted  for.  The  data  in  this  table 
show  that  at  low  angles,  irrespective  of  the  surface  preparation  the 
reflectivities  are  high. 
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1.15  m  laser 


Ge  Diode 


FIGURE  1 


FIGURE  2 
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Introduction 


Successful  predictive  design  of  large  aperture  disc  amplifiers 
requires  a  thorough  knowledge  of  the  operative  physics  in  several 
areas  which  are  not  normally  too  important  in  solid  rod  laser  amplifiers. 
Conversely,  a  disc  amplifier  represents  a  geometry  in  vhich  the  relevant 
physics  can  be  studied  with  a  higher  degree  of  certitude. 

These  effects  of  special  moment  for  disc  amplifiers  can  be  ranked 
in  order  of  importance; 

(1)  free  lasing  and/or  superf luorescence  in  the  discs, 

(2)  since  one  pumps  through  the  faces,  any  active  (i.e.  pumped) 
losses  will  be  present  to  a  much  more  severe  extent  than  in  a  rod 
amplifier, 

(3)  damage  to  the  discs  due  to  the  intense  optical  radiation  of 
the  flashlamps. 

During  the  past  nine  month  period  we  have  explicitly  investigated 
the  first  two  categories  of  effects  and  perforce  the  third.  To  assist 
in  this  effort  we  have  employed  several  powerful  computational  tools 
as  well  as  reasonably  standard  experimental  techniques.  The  Monte-Carlo 
optical  pumping  code  (ZAP)  has  been  used  to  predict  the  energy  storage 
and  a  modified  code  used  to  perform  superfluorescence  calculations, 
quasi- analytic  models  have  been  used  to  predict  free  lasing  thresholds. 

Experimentally,  the  VD  640  glass  laser  was  first  used  for  gain 
measurements  both  with  long  (30  nsec)  and  short  (250  psec)  pulses. 

This  was  not  too  satisfactory  for  several  reasons;  the  spatial  profile 
was  not  too  uniform  and  the  other  demands  for  the  use  of  the  laser 
system  made  it  difficult  to  prosecute  the  experiments  with  as  much 
vigor  as  might  have  been  desired.  The  delivery  of  a  Chromatix  Model 
1000  C  Nd : YAG  laser  in  December  eased  this  restriction  and  additionally 
offered  much  greater  flexibility  in  the  type  of  experiments  which 
could  be  pursued. 
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Free  Lasing  and  Superfluorescence  Results 


Early  experiments  vith  the  French  laser  as  the  pump  laser  showed 
no  gain  anomolles  as  long  as  the  rough  ground  edges  of  the  discs  were 
isolated  from  the  gold  plated  holder  by  a  black  anodized  copper  strip 
(a  first  attempt  to  paint  the  edges  black  was  unsuccessful  since  the 
paint  disintegrated  in  a  few  shots).  The  operation  was  restricted  to 
<  55%  of  the  maximum  bank  energy  by  the  facts  that  we  had  only  15  lamps 
in  these  tests  and  also  that  the  trigger  isolation  capacitors  then 
used  disintegrated  at  18-19  kV. 

In  December  operation  with  the  Chromatix  laser  and  a  full  system 
revealed  that  the  gain  did  not  materially  increase  above  the  level 
previously  reached  as  the  pump  energy  was  increased  to  150  kJ.  This 
was  not  expected  at  that  level  because  of  free  lasing  or  superfluor^ 
escence,  so  a  set  of  experiments  was  prosecuted  to  determine  whether 
or  not  some  active  loss  was  dominating.  This  investigation  is  described 
in  the  next  section  in  more  detail.  It  had  convincingly  negative 
results  which  left  only  free  lasing  and/or  superfluorescence  as  the 
possible  mechanisms. 

Several  alternate  edge  treatments  had  been  under  consideration 
for  some  time  and  a  black  solder  glass  (G9000)  developed  by  Owens- 
Illinois  was  chosen  as  a  more  likely  candidate  than  the  two  ceramic 
glazes  which  were  also  formulated  for  several  reasons: 

(1)  a  lower  firing  temperature  resulting  in  less  chance  of 
having  to  reanneal  the  discs; 

(2)  an  index  of  refraction  more  nearly  matching  the  index  of 
ED-2  (1.7  vs  1.56). 

This  coating  was  applied  to  eight  (8)  discs  and  an  identical 
sequence  was  run  from  25  -  72  7»  of  the  bank  energy  using  the  Chromatix 
laser  to  produce  a  1  kW,  200  nsec  probe  pulse  every  20  jusec  during  the 
pump  pulse.  Figure  1  shows  the  peak  small  signal  gain  coefficient/ 
disc  as  a  function  of  the  edge  treatment.  Two  subsidiary  experiments 
are  included  one,  with  pyrex  glass  shields  in  place  and  one  with  the 
to  lack  copper  strips  removed.  In  the  former  case  the  peak  gains  were 
essentially  the  same  with  and  without  the  shields;  in  the  latter 
case,  fluorescence  probes  revealed  that  even  at  the  lowest  operating 
voltage  at  which  the  lamps  would  trigger  (10  kV) .  the  rolloff  occurred. 
By  extrapolating  the  fluorescence  back  to  a  well  behaved  region  the 
estimate  shown  in  Fig.  1  was  obtained. 

The  analysis  codes  strongly  suggested  that  the  problem  was  relaxa¬ 
tion  oscillations  in  the  discs  rather  than  superfluorescence  and  the 
relative  levels  at  which  oscillation  commenced  with  and  without  the 
black  copper  strips  tended  to  strongly  reinforce  this  conclusion. 

Figures  2  and  3  show  data  which  compares  the  gain  as  a  function  of  time 
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for  the  black  glass  coated  discs  with  the  case  with  copper  strips.  The 
fluctuations  in  the  gain  values  were  marginally  statistically  signi¬ 
ficant  and  much  more  suggestive  of  the  energy  storage  in  an  oscillator 
near  threshold  than  a  system  being  depumped  by  super fluorescence. 

The  results  with  the  black  glass  edge  are  in  excellent  agreement 
with  ZAP  and  the  superfluorescence  code  predicts  a  rolloff  of  only 
~  37.  (which  is  included  in  Figure  1)  due  to  super fluorescence  at  full 
energy.  Figure  4  shows  that  the  full  data  set  with  the  black  glass. 

It  is  also  predicted  that  the  discs  will  not  oscillate  below  e^  gain 
with  this  coating.  There  are  however  two  caveats  at  this  time;  the 
coating  appears  to  be  damaged  to  some  degree  by  ultraviolet  from  the 
flashlamps  and  perfect  adhesion  of  the  coating  to  the  discs  was  not 
found  so  the  full  predicted  gain  may  not  be  feasible.  The  first 
problem  is  being  investigated  now  and  the  second  is  relevant  to 
larger  systems;  this  system  will  operate  with  the  prescat  coating. 

Active  Loss  Investigation 

Earlier  results  on  rod  amplifiers  at  NRL,  AFWL  and  AO  had  suggested 
the  possibility  of  active  losses  playing  some  role  at  high  pump  energy 
densities  and  when  the  gain  rolloff  was  discovered  this  possibility 
was  pursued. 

Operation  with  pyrex  shields  revealed  the  same  rolloff  occurring 
at  essentially  the  same  level  and  ruled  out  ultraviolet  lamp  emission 
below  3000  JL  A  concurrent  experiment  in  which  the  transverse  gain 
through  the  center  of  the  disc  was  measured  showed  the  same  roll  off 
and  argued  against  the  cerium  in  the  ED-2  playing  any  dominant  role 
since  then  it  would  have  been  a  surface  effect  (because  of  the  large 
absorption  cross-section)  and  the  gain  through  the  center  of  the  disc 
would  have  continued  to  rise. 

The  possibility  of  a  broad  band  transient  color  centering  was 
ruled  out  by  spectrally  tuning  the  Chromatix  laser  to  1.12  mp  and 
.946  mp.  In  both  cases  there  was  no  effect  to  a  precision  of  ±.002  cm  . 
This  left  only  the  possibility  of  an  excited  state  absorption  in  the 
neodymium  system  itself.  This  was  checked  by  performing  the  small 
signal  gain  measurement  at  1.052  p,  1.064  p,  1.074  p  and  1.079  p  and 
comparing  the  relative  magnitude  of  the  gain  coefficients  at  the  other 
wavelengths  to  q6  u  an<*  norma^^zin8  the  measurements  to  Duston's 
fluorescence  data’ on  Ed-2.  This  is  shown  in  Figure  5.  In  all  cases 
the  gain  rolled  over  at  ~  507.  pump  energy.  Figure  4  shows  this  data. 

At  the  inversion  level  achieved  (~  .42  J/cm^)  this  data  and  the  highest 
level  data  shows  that  if  there  is  any  excited  state  absorption  it  must: 

(1)  have  a  spectral  shape  very  similar  to  the  fluorescence  curve 
which  would  be  fortuitous,  or 

(2)  have  a  peak  magnitude  <  17./ cm. 
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It  appears  that  an  excited  state  absorption  (if  it  occurs),  is  not  a 
particularly  significant  mechanism  for  disc  amplifiers  since  high  energy 
tensities  compared  to  those  projected  are  required  for  it  to  play  a 
significant  role.  It  should  be  noted  that  there  is  no  basic  inconsis¬ 
tency  between  this  conclusion  and  observations  on  rod  amplifiers  as 
they  were  at  much  higher  stored  energies.  The  superlinearity  of  the 
earlier  NRL  and  AFWL  results  however  suggests  that  they  were  contamin¬ 
ated  by  superfluorescence. 

Damage  to  Disc  Surfaces 

In  conversations  with  LASL  representatives  and  John  Meyers  of 
Owens  Illinois  it  appears  that  LASL  has  seen  massive  surface  damage 
to  discs  and  that  one  possibility  is  that  the  cerium  oxide  powder 
Owens  Illinois  favors  for  final  polishing  forms  an  amalgam  with  the 
glass  with  a  very  high  ultraviolet  absorption  coefficient.  In  the 
course  of  reviewing  why  we  had  not  seen  such  an  effect  we  found  that 

(1)  eleven  of  the  sixteen  NRL  discs  were  repolished  with  Linde  B 
abrasive  after  the  black  paint  disaster, 

(2)  the  remainder  had  to  be  repolished  after  a  subsequent  flash- 
lamp  explosion. 

In  approximately  350  shots  to  date  on  the  disc  amplifier  we  have 
observed  the  following  phenomena  on  a  less  spectacular  but  persistent 
fashion, 

(1)  any  unremoved  polishing  compound  will  be  burned  into  the 
surfaces, 

(2)  we  have  had  a  problem  with  gold  from  the  mounts  flecking 
off  and  being  deposited  on  the  surfaces.  This  seems  to  have  been 
caused  by  improper  surface  preparation  before  plating  and  after  ~  200 
shots  this  effect  has  disappeared, 

(3)  in  several  discs  Internal  bubbles  have  appeared  which  are 
not  oriented  along  the  optic  axis  but  rather  follow  skew  paths  rather 
as  if  they  were  there  initially  and  irradiation  has  Increased  the  size 
to  a  noticable  extent, 

(4)  random  small  damage  spots  apparently  caused  by  dust  particles. 

None  of  these  effects  have  had  any  significant  effect  on  the  laser 
performance  of  the  device  but  they  are  cosmetically  unpleasing. 

Summary  of  Results  to  Date 

The  dominant  problem  in  the  operation  of  the  NRL  system  \as  been 
shown  to  be  pre- lasing  of  the  discs  off  the  disc  edge.  It  appears  that 
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Che  black  glass  coating  Owens  Illinois  has  developed  will  suffice  co 
permit  operation  of  this  device  at  up  to  the  design  levels. 

Projected  Experimental  Program  (6/30/72) 

(1)  Check  geometrical  sealing  of  ZAP  program. 

(2)  Check  gain  to  full  pump  energy. 

(3)  Look  for  excited  state  absorption  at  higher  Inversions. 

(4)  Run  to  full  energy  (400  -  500  J)  with  250  psec  pulses. 

(5)  Measure  pump  uniformity  by  doing  x-y  scan  with  Chromatlx  laser. 
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ABSTRACT 


In  disc  laser  amplifiers,  the  achievable  inversion  (and  thus  the  gain  and 
energy  output)  is  limited  by  fluorescence  amplification  and  parasitic 
oscillation.  In  this  paper,  analytic  results  and  Monte  Carlo  calculations 
of  both  these  effects  are  presented  for  spheres,  circular  discs,  and 
elliptical  discs.  The  effects  of  gain,  fluorescence  line  profile,  refractive 
indices,  and  disc  thickness  on  fluorescence  amplification  are  presented. 

The  loss  rate  is  seen  to  increase  rapidly  above  an  across-diameter  gain  of 
exp  (5)  or  so,  independent  of  thickness.  Parasitic  oscillation  is  analyzed 
as  a  function  of  edge  and  face  index  ratios,  and  of  edge  roughening.  Above 
exp  (3)  gain,  parasitic  suppression  is  difficult,  even  with  a  rough  edge. 
Parasitics  thus  appear  to  present  more  difficulties  than  fluorescence 
amplification.  The  effects  of  both  processes  on  peak  inversion  achieved 
during  a  pumping  pul^  are  also  presented.  In  conclusion,  the  maximum 
practical  size  of  Nd  glass  laser  discs  is  estimated  to  be  about  30  cm 
due  to  the  effects  of  parasitic  oscillation  and  fluorescence  amplification. 
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This  is  a  final  report  on  one  phase  of  the  problem;  work  is  continuing  on  other 
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I.  INTRODUCTION 


Multi-disc  laser  amplifiers  are  presently  of  interest  in  the  production 
of  high-energy  light  pulses  for  the  generation  of  dense,  high-temperature 
plasmas  for  x-ray  generation  and  controlled  thermonuclear  reactions.  Disc 
amplifiers,  as  compared  to  rod  amplifiers,  have  the  advantages  of  more 
uniform  pumping,  easy  heat  removal,  small  size  of  individual  glass  pieces, 
and  less  susceptibility  to  self- trapping.  However,  the  design  of  these 
amplifiers  must  be  carefully  carried  out  if  costly  and  time-consuming  mis¬ 
takes  are  not  to  be  made. 

One  of  the  problems  which  must  be  understood  before  disc  amplifier 

(1-9) 

design  is  undertaken  is  that  of  fluorescence  amplification  (sometimes 

called  amplified  spontaneous  emission,  super  fluorescence,  or  -  incorrectly 
-  superradiance).  Fluorescence  amplification  makes  energy  storage  in 
lasers  more  and  more  difficult  as  the  inversion  level  rises.  Section  II 
of  this  paper  is  a  discussion  of  analytical  and  Monte  Carlo  calculations 
of  the  limitations  due  to  fluorescence  amplification  in  various  geometrical 
shapes  relevant  to  the  disc  laser  problem,  and  of  the  effects  of  varying 
size,  shape,  gain,  refractive  index,  and  line  profile. 

Another  problem  which  arises  in  disc  amplifier  design  is  parasitic 
oscillation.  Such  oscillation  sets  a  sharp  upper  limit  to  the  achievable 
inversion,  and  often  is  a  more  serious  problem  than  fluorescence  ampli¬ 
fication.  In  Section  III  of  this  paper,  parasitic  oscillation  limits  to 
achievable  gain  in  a  single  disc  will  be  calculated  as  a  function  of  the 
reflections  due  to  the  refractive  indices  of  a  disc  and  its  surroundings, 
and  as  a  function  of  the  edge  roughness  of  the  disc. 
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Both  fluorescence  amplification  and  parasitic  oscillation  act  to  limit 
the  maximum  practical  size  of  the  discs  in  a  disc  amplifier,  since  both  pro¬ 
cesses  get  worse  as  a  disc  becomes  larger  (assuming  constant  inversion  density). 
Accurate  analysis  of  both  effects  is  therefore  required  to  determine  the  exact 
limitations  on  size  and  gain. 

II-  FLUORESCENCE  amplification 

Many  materials,  wnen  properly  pumped,  exhibit,  fluorescence  (the  sponta¬ 
neous  emission  of  photons).  The  fluorescence  arises  when  the  upper  level  of 
the  fluorescent  transition  is  populated  by  the  pumping  process,  and 
radiative  transitions  to  the  lower  level  take  placo.  Once  the  pumping  is 
turned  off  the  fluorescence  usually  decays  exponentially  with  a  time 
constant  called  the  fluorescent  lifetime  T.  An  exponential  decay  implies 
that  the  process  is  linear  -  that  is,  that  the  fluorescence  is  linearly 
proportional  to  the  population  in  the  upper  level. 

Suppose,  however,  that  the  pumping  was  so  intense  that  an  inversion 
exists;  the  population  in  the  upper  level  is  greater  than  the  population 
in  the  lower  level.  The  pumped  material  then  exhibits  optical  gain  at  the 
transition  frequency,  and  may  be  used  to  amplify  light  in  a  laser  amplifier 
or  oscillator.  However,  the  fluorescence  continues  under  conditions  where 
gain  exists;  in  fact,  the  processes  of  spontaneous  emission  (fluorescence) 
and  stimulated  emission  (gain)  are  inextricably  interconnected.  This  means 
that  fluorescent  light,  once  emitted  within  the  volume  of  the  material, 
will  be  amplified  before  it  reaches  the  edges.  This  increases  the  effective 
fluorescence  loss  rate,  which  is  the  same  as  decreasing  the  instananeous 
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fluorescent  lifetime.  Because  the  gain,  and  thus  the  added  loss,  now  depends 
on  the  inversion  level,  the  fluorescent  decay  process  is  no  longer  lineal. 
Instead,  the  decay  rate  becomes  faster  as  the  inversion  level  increases. 

This  fluorescence  amplification  problem  is  especially  severe  in  large  laser 
systems,  since  the  long  path  lengths  available  lead  to  large  total  gain  and 
thus  large  amounts  of  fluorescence  amplification. 

So  far  we  have  considered  what  happens  when  the  pumping  is  over. 

However,  the  fluorescence  and  gain  processes  also  exist  during  the  pump 
pulse.  If  the  fluorescence  amplification  is  negligible,  so  that  the  system 
is  linear,  then  doubling  the  pumping  amplitude  will  double  the  inversion. 

In  general,  the  peak  inversion  will  be  linearly  proportional  to  the  pump 
pulse  amplitude.  This  proportionality  will  be  destroyed  if  gain  exists  in 
the  material.  In  the  presence  of  fluorescence  amplification,  equal  incre¬ 
ments  of  pumping  will  yield  smaller  and  smaller  inversion  (and  gain) 
increases,  because  the  decay  rate  per  unit  of  inversion  will  increase  with 
the  inversion.  The  details  of  pumping  in  the  presence  of  fluorescence 
amplification  are  studied  in  Section  IV. 

In  this  section  we  will  consider  only  single-pass  fluorescence 
amplification.  That  is,  the  edges  of  the  discs  will  be  assumed  to  be 
totally  absorbing,  so  that  once  the  fluorescent  light  hits  an  edge  it  is 
not  reflected  back  into  the  disc  again.  This  is  the  best  possible  edge 
condition  for  an  amplifier,  since  if  the  disc  edge  is  not  perfectly  black 
the  light  will  be  reflected  back  into  the  disc  and  excite  further  stimulated 
emission,  thus  increasing  the  fluorescence  loss  rate.  We  are  therefore 
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calculating  an  upper  bound  on  gain  due  to  the  effects  of  fluorescence 
amplification,  since  any  actual  disc  will  return  some  reflected  light,  thus 
increasing  the  loss  rate  and  decreasing  the  gain.  We  here  limit  ourselves 
to  the  black-edge  case  because  once  light  is  reflected  from  the  edges 
oscillation  will  occur  at  some  gain  level  and  a  different  method  of  analysis 
must  be  used  to  find  this  level  (see  Section  III). 

A.  Fluorescence  Amplification  in  a  Sphere 

We  first  calculate  the  stimulated  emission  in  a  sphere,  not  because 
sphere  lasers^  ^  are  of  interest  but  because  an  exact  analytical  calculation 
is  possible  in  this  case.  In  addition,  the  sphere  forms  a  simple  system  in 
which  to  demonstrate  the  effects  of  line  shape. 

1.  Single-Wavelength  Calculation 

Consider  a  laser  material  in  the  shape  of  a  sphere.  We  will  assume  that 

a  spatially  uniform  inversion  exists  throughout  the  sphere,  and  make  a 

calculation  of  the  instantaneous  rate  of  fluorescence  amplification.  The 

inversion  causes  the  material  to  have  a  gain  coefficient  a,  such  that  the 

intensity  of  a  pencil  of  light  rays  increases  as  P  <=  Pq  exp  (a-l),  where  l 

is  the  path  length  along  the  ray  direction.  Suppose  that  spontaneous  emission 

releases  a  power  of  I  per  unit  volume.  If  the  sphere  has  a  diameter  D,  the 

3 

total  spontaneous  emission  is  I  it  D  /6.  We  will  calculate  the  amount  by 
which  this  spontaneous  emission  is  amplified. 

Introduce  a  system  of  spherical  coordinates  centered  on  the  surface  of 
the  sphere  (Fig.  1).  The  angle  0  is  measured  from  the  sphere  diameter, 
and  0  gives  the  rotation  around  the  diameter.  The  distance  r  is  measured 
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Fig.  1  -  Coordinate  system  for  the  analysis  of  fluorescence 
amplification  In  a  sphere  of  diameter  D.  The  origin  Is  at  the 
sphere's  south  pole. 


along  the  direction  defined  by  9  and  0.  The  element  of  volume  In  these 
o 

coordinates  is  r  sin0  dr  d0  d0,  and  the  fluorescence  from  this  volume  is 
I  r2  sin0  dr  d9  d0.  We  assume  that  this  fluorescence  is  emitted  isotropically. 
An  element  of  sphere  surface  of  area  dS  at  the  origin  of  coordinates  will 
intercept  a  fraction  of  the  fluorescence  equal  to  the  projected  area  of  dS  in 
the  direction  of  the  volume  element  (dS  cosq,)  divided  by  the  area  of  a  sphere 
centered  on  the  volume  element  and  passing  through  the  origin  (4  n  r2). 

However,  the  fluorescence  will  be  amplified  by  a  factor  exp  (a?)  before  it 
reaches  dS.  Thus  dS  will  intercept  a  total  radiation  (spontaneous  plus 
stimulated)  of 


exp  (or)  cas<t>  sin 0  dr  d8  d0  dS. 


Now  integrate  this  expression  over  the  volume  of  the  sphere  to  find  the 
total  flux  on  the  area  dS: 

2  it  n/2  D cas<J> 


dP 


TT  /d9/d*  / 


dr  e0*  cas<t>  sin  0 


o  o  o 


n/2  Datf0 


Vs  /-*/ 


dr  e0*  «W0  sin0 


o  o 
n/2 


=  j  d0  ^ea  j  OTT0  Fin0  . 


F-10 


With  the  change  of  variable  u 


=  CVS $,  the  flux  becomes 


1 

o 


I  ds 

2a 


By  symmetry,  all  elements  of  the  sphere  surface  receive  the  same  flux,  so  the 
total  emission  hitting  the  sphere  surface  is 


To  find  the  amount  M  by  which  the  original  spontaneous  emission  has  been 

3 

multiplied,  divide  by  the  total  spontaneous  emission  I  n  D  6  and  find 


where  we  have  introduced  the  variable  0  -  oP.  Note  that  0  is  the  log  of  the 
gain  straight  across  the  sphere.  The  ratio  of  stimulated  to  spontaneous 

emission  is 


A  = 


3_ 

2p 
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Our  calculation  has  been  for  a  single  value  of  a«  However,  in  actual 
laser  materials  the  fluorescence  and  gain  are  distributed  in  wavelength 
according  to  some  line  profile,  with  values  from  zero  to  the  peak  value  being 
represented.  The  effective  values  of  the  spontaneous  emission  multiplier  M 
and  the  stimulated-spontaneous  ratio  A  are  therefore  found  by  averaging  the 
single  wavelength  values  over  the  line.  From  the  viewpoint  of  fluorescence 
amplification,  the  worst  possible  case  would  be  a  flat-topped  or  rectangular 
line  profile,  since  in  that  case  the  line  average  is  equal  to  the  peak  value. 
For  more  realistic  line  profiles,  in  which  values  less  than  the  peak  are 
represented,  the  line  average  is  less  than  the  peak  a  by  an  anount  which 
depends  on, the  line  shape. 

It  should  be  noted  that  M  =  1  +  A  for  any  line  shape,  since  by  definition 
M  =  (J  +  K)/J  =  1  +  J/K  and  A  =  J/K,  where  J  is  the  spontaneous  radiation 
and  K  is  the  stimulated  radiation. 

2.  Lorentzian  Line  Profile 

Let  us ,  for  example,  calculate  the  values  of  M  and  A  for  a  Lorentzian 
line  profile.  The  gain  (and  the  fluorescence,  which  has  the  same  line 
shape)  have  the  form 


a  = 


l  + 


-  V 

w 


where  is  the  gain  at  the  peak  of  the  line,  \q  is  the  center  wavelength  and 
W  is  the  full  width  at  half  maximum. 
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Since  all  that  matters  when  averaging  over  a  line  is  the  fraction  of  the 
line  representing  different  values  of  a>  we  need  to  find  a  function  which  gives 
this  amplitude  weighting.  For  any  symmetric  shape  of  a,  we  have 


M  = 


l 


+  00 

a  M  dX 


/ 


+® 
a  dX 


f  a  M  da  da 

0 

/a 

e 


Introducing  the  variable  x  *  aJ Op>  we  have 

,1 

M  =  J  M  (o^)  (^(x)  dx 

o 

where  the  Lcrentzian  weighting  function  (^(x)  is  given  by 


a 


dX 


plu>  ■ 

'  Vii 


/^d“ 


.  >/x(l-x) 


f  dx 
J  Vx(l-x) 
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. jig.  ...  . . 


- '  *■*"-  ■ 


Since  the  integral  in  the  denominator  is  the  beta  function  =  n, 

we  have 


PL(x)  - 


JT  VX(l_X) 

The  line-averaged  value  of  the  sli  iulated/ spontaneous  ratio  is  given  by 


7  A  (v) 


PL(x)  dx 


1 

3_  f 

2ePx  ( 

\  JL\+_L.  _1 

i 

2n  J  dx 

o 

S2x2  \ 

; "  *)  p3>c3 "  p* 

. 

Vx(i-x) 

where  3  =  cx^D.  Unfortunately,  the  various  terms  of  the  integrand  lead  to 
(cancelling)  infinities.  We  therefore  calculate  AL  by  a  less  direct  method. 

Firstly,  we  determine  the  power  series  expansion  of  A  by  expanding  the 
exponential  and  combining  terms.  We  find 

3A  .V  .hL.^L 


A  8  10  48  280 


3P 


k. 


n=  1 


(nf3) (nfl) ! 


where  »  aP. 
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Then  the  line  average  is 


pV‘* 


(l-x)2  (irf3)  (n+1) ! 


We  reverse  the  order  of  summation  and  integration,  and  use  the  fact  that  the 
integral  over  x  is  a  beta  function  to  find 


A  .  2B+  3fi^+  5fl?  +  ^_  + 
\  16  80  768  1024 


_ 

(n+3) (n+1) ! 


n=l 


n 

n 

r=l 


Since  the  series  must  be  summed  numerically  to  find  the  corresponding  to 
some  specific  p,  the  error  due  to  truncating  at  any  n  must  be  known.  Let 
the  relative  error  due  to  ignoring  terms  beyond  m  be  so  that. 

v-fe  ^")  K) 

where 

_  3  /  2r-l\ 

an  “  (n+3)  (nf  1)  !  Jj=1  \  2r  /  * 
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Then  we  have 


As  m  becomes  large,  the  ratio  of  successive  coefficients  of  (3  in  the  expansion 
of  A,  becomes 


anfl  =  (rri-3)  (2nf 1)  1 

a  (nf4)(nf2)(2nf2)  7 

r.  n  +  — 


We  therefore  have  an  upper  limit  to  the  truncated  part  of  the  series: 


V  fin  V  a"*1  /  B  \n  _  V 
^  anP  <  E n  amflP  (mf  9/ 2  ) 

n=mfl  n=0  \  /  1  - 


.r1 


mf  9/  2 


Therefore 


€  <-/ 
m 


nrt- 


H)e 


m 


where  Em  =  E  arf' 

n=l 


3.  Gaussian  Line  Profile 


We  may  find  the  effect  of  a  Gaussian  line  profile  having  the  form 


a  = 


cl 
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by  the  same  method  used  for  the  Lorentzian  line.  In  this  case,  the  amplitude 
weighting  function  is 


pg(x)  -  ~r~=r 

«/-  n  x 

and  term-by-term  integration  of  the  series  expansion  of  A  yields 


+  JLi 

30 


96 


280/5 


-f  _ & _ _ 

^  (r*t3)(nfl)I 
n=l 

In  this  case,  as  in  the  Lorentzian  one,  the  ratio  of  successive  coefficients 
approaches  n  +  7/2  as  n  -  ®,  and  so  the  truncation  error  limit  given  for  the 
Lorentzian  line  profile  also  applies  to  the  series  expansion  of  A  for  a 
Gaussian  line  profile. 

4.  Summary  of  Analytic  Sphere  Results 

Table  I  gives  numerical  values  of  the  ratio  of  stimulated  to  spontaneous 
radiation  A  as  a  function  of  (3  (the  log  of  the  across-sphere  gain  at  the 
line  peak)  for  flat-topped,  Gaussian,  and  Lorentzian  line  profiles  of  the  gain 
coefficient.  For  small  values  of  3,  we  see  that  A  is  linearly  dependent  on 
3  with  a  coefficient  which  depends  on  the  line  shape.  The  rise  of  A  with  3 
becomes  faster  as  3  approaches  unity,  and  is  very  fast  for  values  of  3  much 
greater  than  one.  We  shall  see  later  that  this  behavior  is  typical  of  any 
shape  of  laser  material.  The  fact  that  decay  rates  rise  sharply  for  g  >  1 
means  that  pumping  becomes  more  and  more  difficult  in  this  region. 
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TABLE  I 


Ratio  of  stimulated  to  spontaneous  radiation  A  in  a  sphere  as  a  function  of 
the  log  of  the  across- sphere  gain  at  the  line  peak  j3  for  flat-topped,  Gaussian, 
and  Lorentzian  line  profiles. 


a 

^LAT 

agauss 

^LORENT 

.01 

.001760 

.002657 

.001879 

.02 

.007540 

.005326 

.003765 

.05 

.01900 

.01340 

.009470 

.1 

.03852 

.02710 

.01913 

.2 

.07917 

.05543 

.03905 

.5 

.2153 

.1484 

.1040 

1. 

.5000 

.3352 

.2326 

2. 

1.396 

.8790 

.5980 

5. 

12.97 

6.49 

4.153 

10. 

593.6 

206.8 

123.4 

20. 

3.457  x  106 

8.097  x  105 

4.670  x  105 

50. 

6.097  x  1018 

8.778  x  1015 

4.992  x  1017 

100. 

39 

7.984  x  10 

QO 

8.055  x  10 

4.562  x  1Qj8 

5.  Monte  Carlo  Calculation 

The  amount  of  fluorescence  amplification  in  a  sphere  was  also  calculated 
by  means  of  a  Monte  Carlo  optical  power  flow  program.  Since  analytical 
results  are  available  in  the  case  of  a  sphere,  it  may  seem  pointless  to 
duplicate  the  results  with  an  approximate  method.  However,  the  Monte  Carlo 
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program  can  be  used  with  geometric  shapes  of  arbitrary  complexity.  Thus  it 
is  much  more  general  in  its  application  than  analytic  methods,  but  of  course 
any  such  program  should  be  checked  in  a  geometry  where  the  results  are  known. 
This  is  the  reason  for  its  use  in  the  case  of  a  sphere. 

The  Monte  Carlo  program  is  called  UAP.  It  traces  optical  power  flow  in 
systems  which  are  made  up  of  segments  bounded  by  planes,  spheres,  cylinders . 
cones,  conic  surfaces,  and  so  forth,  or  in  combinations  of  such  segments.  The 
material  in  the  segments  is  specified  by  its  refractive  index  and  absorption 
(or  amplification)  coefficient.  The  system  surfaces  may  have  constant  or 
angle-varying  reflection,  both  specular  and  scattering.  Each  optical  ray 
carries  with  it  up  to  two  hundred  wavelength  intervals.  The  bulk  and 
surface  properties  of  the  system  may  vary  with  wavelength,  and  the  optical 
power  in  each  wavelength  interval  will  change  according  to  the  system 
properties  at  its  particular  wavelength.  Rays  are  started  at  random  from 
specified  sources  (surface  or  volume)  with  given  initial  wavelength 
distributions  of  power,  and  traced  through  the  system  to  find  the  deposition 
or  extraction  of  power  in  all  the  system  elements.  A  ray  cut-off  scheme  is 
used  which  assures  that  ray  truncation  error  is  small  without  wasting  time 
following  weak  rays.  ZAP  was  designed  for  accurate  modelling  of  laser  pump 
cavities  (a  job  it  does  well),  but  it  is  also  well  suited  to  the  calculation 
of  fluorescence  amplification  in  arbitrary  geometries. 

The  sphere  geometry  was  specified  to  ZAP,  and  runs  of  2500  rays  were  made 
for  various  values  of  (5  to  achieve  an  accuracy  of  two  or  three  percent. 

The  rays  were  started  at  random  throughout  the  sphere,  to  correspond  to 
uniform  inversion.  All  rays  hitting  the  sphere  surface  were  totally  absorbed. 
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Runs  were  made  for  a  flat-topped  line  (one  wavelength),  for  a  line  with  a 
Gaussian  profile,  and  for  a  line  with  a  Lorentzian  profile.  The  Gaussian  and 
Lorentzian  lines  ware  modelled  by  dividing  a  line  profile  into  equal  intervals 
of  wavelength.  The  gain  coefficient  and  starting  power  in  each  wavelength 
interval  were  then  taken  proportional  to  the  value  of  the  gain  coefficient  at 
the  center  of  the  interval.  It  is  necessary  to  go  a  large  distance  out  on 
the  tail  of  the  Lorentzian  line  in  order  to  get  a  proper  weighting  of  low 
values  of  the  gain  (recall  that  the  Lorentzian  amplitude  weighting  function 
approaches  infinity  as  the  amplitude  approaches  zero).  The  final  result, 
with  the  above  number  and  spacing  of  wavelengths,  is  about  1%  high  for  a 
sphere  when  the  calculation  is  exact.  To  this  systematic  error,  the 
Monte  Carlo  method  adds  the  previously  mentioned  random  error  of  2  or  3%. 

A  comparison  of  the  analytic  and  Monte  Carlo  results  for  a  sphere  is 
shown  in  Fig.  2.  The  agreement  is  quite  satisfactory,  since  the  exact 
result  lies  within  the  Monte  Carlo  error  limits  in  all  cases. 

B.  Circular  Disc 

We  will  now  consider  the  case  of  a  circular  disc,  which  is  of  considerably 
greater  practical  interest.  We  will  first  perform  an  approximate  analytic 
calculation  which  is  exact  for  infinity  thin  discs  at  low  gain,  and  then 
find  the  results  for  thick  discs  and/or  high  gain  by  the  Monte  Carlo  method. 

1.  Analytic  Approximation:  Thin  Disc,  Low  Gain 

Consider  the  spontaneous  fluorescence  emitted  from  a  small  element  of 
volume  in  a  circular  disc  with  a  totally  absorbing  edge.  Some  of  the  radia¬ 
tion  will  strike  the  edge  directly,  while  the  re3t  will  hit  the  faces.  Both 
components  will  be  amplified  as  they  travel.  The  radiation  which  goes  to 
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the  faces  can  be  further  separated  into  two  parts.  Some  will  strike  at  less 
than  the  critical  angle  for  total  internal  reflection,  in  which  case  most 
will  be  transmitted  and  a  small  fraction  (  rv  5%)  reflected  back  into  the  disc. 
Some  will  strike  at  more  than  the  critical  angle,  and  will  be  totally 

reflected  back  into  the  disc.  This  totally  reflected  component  will  continue 

to  bounce  back  and  forth  between  the  faces  until  it  finally  strikes  the  black 
edge.  The  path  length  for  this  process  is  the  same  as  if  the  radiation  had 
gone  straight  to  the  continuation  of  the  edge  beyond  the  faces  (Fig.  3). 

Thus  there  are  two  parts  of  the  circular  disc  gain  calculation:  the 
radiation  below  the  critical  angle,  and  that  above  the  critical  angle, 
a.  Radiation  below  the  critical  angle 

Consider  a  circular  disc  of  diameter  D  and  thickness  L.  A  volume  element 
at  a  distance  r  from  the  disc  center  emits  a  pencil  of  radiation  at  an  angle 
0  from  the  normal  to  the  disc  (Fig.  4).  The  distance  from  the  point  of 

emission  to  the  edge  is  l;  the  ray  travels  a  distance  l  esc  0  before  it  hits 

the  edge.  Note  that  this  result  is  independent  of  the  position  along  an 
axis  normal  to  the  disc  faces  (the  z-position)  of  the  emitter.  The  distance 
t  depends  on  the  angle  8  from  a  disc  radius  through  the  emitting  element  to 
the  emitted  pencil.  From  Fig.  5  we  see  that 


l  = 


(r  sin  9) 2 


r  avQ  . 


The  radiation  is  multiplied  by  a  factor  of  exp  (a  l  esc  0)  as  it  travels  to 
the  edge;  we  recall  that  a  is  the  gain  coefficient  in  the  material. 
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Fig.  5  -  Geometry  for  the  calculation  of  the  distance  from 
the  point  of  emission  to  the  disc  edge.  The  formula  is  the 
same  for  9  <  jt/2  and  9  >  n/2,  as  shown. 


We  must  Integrate  over  the  disc  volume  (i.e.,  over  r)  and  over  solid 
angle.  The  fractional  disc  volume  at  radius  r  is  8r  d’2  dr,  and  the  fraction 
of  solid  angle  at  0,  0  is  (4n)  1  sin  0  d0  d0.  The  angle  0  varies  from  0 

c 

to  n-0c,  where  @c  is  the  critical  angle,  but  by  symmetry  we  can  take  twice 
the  integral  from  9c  to  n/2.  Likewise  we  may  take  twice  the  integral  from 

0  60  I!  ter  The  above-critical- angle  integral  for  the  ratio  M  of  total  to 
spontaneous  emission  is  then 


D/2  n/2  n 

■  — —  £  fir  J'wJ' ^9  r  sin  0 


exp 


0  0  0 
c 


<*(Vy  -  r2  sin2  0  -  r 


sin  0 


This  is  somewhat  clarified  by  use  of  the  dimensionless  variable  u  -  2  r/D, 
whence 


n/2  n 

^0  u  sin  0  exp 

I  0 
c 

op(y  1  -  U2  sin2e  -  U  awe) 

2  sin  0 

l0W  8Sin  condition  aP  «  1,  we  can  expand  the  exponential  in  its 
power  series  and  produce  a  simpler  integral: 

1  n/2  n 


mi  *  i  yv d0/ u 


sin  0 


0  0  0 
c 


1  + 


aP 


(VT7“ 


2  2 

sin  0  -  flare 


. 


2  sin  0 


The  first  term  in  the  series  yields  the  simple  result  but  the  second  term 

is  «.ore  difficult.  Sin  0  cancels  in  numerator  and  denominator,  so  the  0  integral 
is  Just  (n/2  •  0c),  leaving 
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M, 


cos  8 


+  f  (f  -  9c)/du/d'’  u  (V1  -  “ 

'  *  a  r\ 


2  2 

sin  0  -  u  COS  0 


)• 


In  the  0  integral,  the  COS  0  portion  integrates  to  zero,  leaving  twice  the 
complete  elliptic  integral  of  the  second  kind 

n 

2E  =  f  ^  ^1  -  sin^  0 


so  that 


which  by  the  change  of  variable  q  ■  u  becomes 


1 

M1  “  m  0c  +  f  (f  "  9c)  /dq  E(q) 


thus  yielding  the  final  result 


«l’“  Sc  +  f  (f  •  8c)(%+0)  1 


where  G  ■  .9159656  •••  is  Catalan's  constant. 
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£.  Radiation  above  the  critical  angle 


Let  us  follow  the  radiation  emitted  at  above  the  critical  angle  for 
total  internal  reflection  by  a  volume  element  at  distance  b  from  the  disc  face 
(Fig.  6).  Starting  at  angle  0  from  the  normal  to  the  face,  the  ray  travels 
a  distance  b  sec  0  before  it  hits  the  face.  The  relative  strength  is  now 
exp  (a  b  sec  0).  A  fraction  T  is  transmitted  through  the  face  and  (for  our 
purposes)  lost,  while  (1-T)  is  reflected  and  crosses  the  disc.  The  path  length 
to  the  opposite  face  is  L  sec  0  (recall  that  L  is  the  disc  thickness),  and 
so  the  reflected  ray  is  amplified  by  exp  (a  L  sec  0).  It  then  again  loses 
a  fraction  T  to  the  outside,  and  (1-T)  recrosses  the  disc,  is  amplified,  and 
is  partially  transmitted  and  partially  reflected.  The  bouncing  process 
continues  until  the  ray  hits  an  absorbing  edge.  Let  the  relative  strength 
at  the  first  encounter  of  the  face  be  P  =  exp  (a  b  sec  0),  and  define  the 
gain  on  a  pass  all  the  way  across  the  disc  to  be  K  =  exp  (a  L  sec  0).  Then 
the  ray  strength  is  as  follows  (Fig.  7): 


1 

P 

P(l-T) 

P(1-T)K 

P(1-T)2K 

2  2 
P(l-TrK 

3  2 
P(l-T)  K 

3  3 
P(1-T)JK 


at  the  start 

at  the  initial  face  intercept 
after  t.he  first  reflection 
after  crossing  to  the  opposite  face 
after  the  second  reflection 
after  the  second  full  crossing 
after  the  third  reflection 
after  the  third  full  crossing 


and  so  forth. 

+  P2U-T)V  + 


The  total  amount  radiated  is  then  PT  +  PT(1-T)K  +  PT(1-T)2K2 
...  which  sums  to  PT/[1  -  (l-T)K]  if  we  take  an  infinite  number 
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of  terms.  Of  course,  the  actual  ray  will  terminate  after  a  finite  number  of 
bounces,  but  the  error  in  the  infinite-bounce  approximation  is  small  as  long 
as  the  gain  across  the  disc  times  the  reflectivity  is  small,  in  which  case 
the  ray  strength  is  reduced  to  a  negligible  value  in  just  a  few  bounces.  This 
smalx-gain  case  is  exactly  the  one  we  are  calculating,  so  we  will  use  the 
inf ini te- bounce  approximation  above,  which  overestimates  the  actual  result. 

We  must  integrate  over  angle  and  position  to  find  the  overall  average 
gain.  The  fractional  solid  angle  at  angle  0  is  sin  0  d0/2,  and  the  fractional 
volume  is  db/L  (note  that  the  r  and  0  integrals  have  vanished  because  there  is 
no  r  or  0  dependence  in  the  infinite-bounce  approximation).  Then  the  above- 
critical-angle  radiation  contributes  approximately 


M 


2 


0  0 


T  exp (g  b  sec  0)  sin  tb 
1  -  (1-T)  exp (a  L  sec  0) 


where  we  have  (by  symmetry)  taken  twice  the  0  integral  from  0  to  0  rather 

c 

than  the  integral  from  0  to  0  plus  the  integral  from  it  -  0  to  it.  For 

c  c 

«  1  we  expand  the  exponentials  in  the  integrand  and  divide  to  get 


9c  L 

/d0  Jdb  sin  0  ^  1  +  ^ -S-L 

J  j  V  T  COS  *0 

0  0  ' 


+ 
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which  Integrates  to 


While  performing  this  calculation  we  have  made  another  approximation:  we  have 
assumed  that  the  transmission  T  is  independent  of  angle,  and  have  ignored  the 
fact  that  T  depends  on  polarization.  In  actual  fact,  T  decreases  as  $ 
approaches  the  critical  angle,  and  more  radiation  returns  to  the  disc  to  be 
amplified.  If  we  use  the  normal-incidence  value  for  T,  this  approximation 
underestimates  the  actual  result,  and  (hopefully)  tends  to  offset  the  over¬ 
estimate  due  to  the  infinite-bounce  approximation.  In  any  case,  both 
approximations  get  better  as  the  disc  gets  thinner,  since  then  even  a  ray 
starting  near  the  edge  takes  many  bounces  to  get  to  that  edge,  and  the 
across-disc  gain  is  low  so  even  if  a  ray  takes  more  bounces  to  exit  than 

we  calculate  (due  to  greater  reflectivity)  the  added  radiation  is  small. 

Y.  Total  Result 

Adding  together  the  results  above  and  below  the  critical  angle,  we  find 
that  for  low  gain  (aD  «  1  and  at  «  1)  and  thin  discs  we  have 

„,1  +  0D(l.^(i  +  G)+aL(l. 

This  result  is  for  one  value  of  a,  and  thus  corresponds  to  a  flat  line 
profile.  Values  of  m  corresponding  to  other  line  shapes  may  be  calculated 
as  in  the  sphere  case. 
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With  a  Lorentzian  line  profile,  the  amplitude  weighting  function  Is 

We  can  write  M  In  the  form  M  =  1  +  Ba(X),  where 


^(x)  *■  £ic/k(1-x)  J  “1. 

B  Is  Independent  of  X.  Then 


K,  a/to  1  +  B«* 

0  xjx(l-x) 


i  +  Sb  s 

it  2  * 


where  a  Is  now  the  gain  coefficient  at  the  line  peak. 
We  thus  have 


H,  s:  1  +  £  B, 
u  L 


or 


"Lal  +  “° 

With  a  Gaussian  line  profile,  p^,  =  <^r  and  we  get 


M_  —  f  dx  1  +  B  »  x 
G  J  /-it  Bn  x 
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X 


1 

*  l  +  f  dx 

>  { 


J-fa  x 


-  I  +  1=fl- 


We  thus  have 


mg 


B, 


or 


Since  the  low-gain  result  In  any  geometry  will  always  have  the  form 
M-l  +  CaorA^Ca  (where  C  Is  a  constant  depending  on  the  geometry);  the 
.  flat-topped;  Gaussian,  and  Lorentzlan  lines  will  always  have  the  low-gain  ratios 
2 : «/2 : 1  exhibited  for  the  sphere  and  cylinder. 

2.  Monte  Carlo  Calculation 

ZAP  was  used  to  find  the  stlmulated-spontaneous  ratio  A  In  circular 
discs.  Again  the  Inversion  (and  thus  the  gain  coefficient)  was  assumed 
spatially  uniform.  The  reflection  at  the  face  was  taken  as  the  average  over 
polarization.  In  Fig.  8,  we  see  the  results  for  flat,  Gaussian,  and 
Lorentzlan  line  profiles  In  discs  whose  thickness  is  20%  of  their  diameter. 
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Fig.  8  -  Monte  Carlo  values  of  the  ratio  A  of  stimulated  to 
spontaneous  loss  as  a  function  of  the  product  g  of  the  gain 
coefficient  and  diameter  for  a  circular  disc  of  thickness  equal 
to  20%  of  its  diameter.  The  ratio  of  the  refractive  index 
inside  the  disc  to  that  outside  the  faces  is  1.56.  Curves  are 
shown  for  three  different  line  profiles  of  the  gain  coefficient. 
The  low-gain,  thin-diac  approximation  is  shown  for  comparison. 
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The  refractive  Index  Inside  the  discs  was  1.56  times  the  index  outside  and  the 
circular  edge  of  the  disc  was  totally  absorbing.  The  A  ratios  are  shown  as  a 
function  of  p  ■  a  D,  where  a  is  the  gain  coefficient  at  the  line  peak  and  D 
is  the  disc  diameter.  Also  shown  is  the  small-gain,  thin-disc  approximation. 
The  approximation  is  seen  to  be  good  below  p  a  .1,  even  for  discs  this  thick. 

The  A  curves  have  the  same  general  shape  as  the  curves  for  a  sphere, 
except  for  a  faster  rise  when  A  is  greater  than  one  (this  faster  rise  is  due 
to  the  longer  path  available  in  the  cylinder  due  to  total  internal  reflection 
from  the  faces). 

The  aswunt  of  fluorescence  amplification  in  a  disc  depends  on  the  disc 
thickness  and  on  the  ratio  of  refractive  indices  inside  and  outside  the 
sphere,  as  well  as  on  the  gain  coefficient.  As  the  cylinder  index  approaches 
the  index  of  its  surroundings,  less  and  less  of  the  light  inside  it  under¬ 
goes  total  internal  reflection,  and  A  becomes  lower.  For  exact  index  matching, 
no  light  returns  into  the  disc  and  A  has  its  minimum  value.  The  variation  of 
A  with  index  ratio  is  shown  in  Fig.  9  for  various  values  of  p.  The  low- 
gain,  thin-disc  approximation  is  also  shown  in  Fig.  9,  and  is  seen  to  match  the 
Monte  Carlo  results  well  except  at  low  index  ratios,  where  the  angle  for 
total  internal  reflection  increases  and  the  infinite- bounce  approximation 
breaks  down. 

The  Monte  Carlo  results  for  the  variation  of  A  with  relative  disc 
thickness  are  shown  in  Fig.  10.  The  values  for  p  ■  .1  are  compared  to  the 
low-gain,  thin-disc  approximation.  The  approximation  is  seen  to  be  good  to 
at  least  80%  disc  thickness.  The  most  interesting  feature  of  theoe  results 
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Fig,  9  -  The  effect  of  inside-outside  refractive  index  ratio  on 
the  fluorescence  amplification  in  a  20%  thickness  circular  disc. 
The  increase  with  index  ratio  is  due  to  the  greater  fraction  of 
radiation  trapped  by  total  internal  reflection  at  higher  ratio. 
The  low  gain,  thin-disc  approximation  is  compared  to  the 
P  *  .1  result,  which  it  fits  exactly  except  at  low  ratios. 
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FRACTIONAL  THICKNESS 


Fig.  10  -  The  fluorescence  amplification  In  a  disc  of  1.56  inside- 
outside  refractive  index  ratio  as  a  function  of  the  ratio  of  the 
disc  thickness  to  its  diameter.  The  increase  is  slow  because  the 
principal  contribution  is  from  rays  undergoing  total  internal 
reflection,  and  this  contribution  is  independent  of  thickness.  The 
low-gain,  thin-disc  approximation  is  compared  with  the  p  -  .1  result. 
Note  broken  vertical  scale. 
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is  the  very  slow  increase  of  fluorescence  amplification  with  disc  thickness. 
This  slow  increase  arises  because  much  of  the  fluorescence  amplification  in  a 
disc  is  due  to  rays  which  are  totally  reflected  at  the  faces  of  the  disc,  and 
the  contribution  of  such  rays  is  independent  of  disc  thickness.  Since  useful 
laser  gain  is  roughly  proportional  to  disc  thickness,  discs  should  be  as 
thick  as  possible  within  the  constraints  of  pumping  uniformity,  thermal 
distortion,  self-trapping,  and  fabrication  difficulty.  There  is  also  an 
obvious  advantage  in  trading  a  thick,  weakly-doped  disc  for  a  thin,  highly- 
doped  disc,  since  fluorescent  lifetime  usually  decreases  with  doping. 

C.  Elliptical  Disc 

Analytical  computation  is  difficult  in  the  case  of  an  elliptical  disc, 
and  unnecessary  in  light  of  the  good  fit  of  the  Monte  Carlo  results  to  the 
sphere  and  circular  disc  cases,  Therefore  only  Monte  Carlo  calculations 
were  made  in  this  case.  The  computations  were  for  discs  whose  major  and  minor 
axes  were  in  the  ratio  1:0.5,  which  is  typical  of  laser  glass  tilted  at 
Brewster's  angle  in  air  if  the  beam  is  round.  The  thickness  was  1/7  of  the 
major  axis,  since  this  corresponds  to  an  actual  disc  in  use  in  a  multidisc 
laser  amplifier  in  oui  laboratory.  Once  more,  the  inversion  was  uniform, 
the  inside-outside  refractive  index  ratio  was  1.56,  and  the  curved  edge  was 
totally  absorbing. 

The  ratio  of  stimulated  to  spontaneous  radiation,  A,  is  shown  in  Fig.  11 
as  a  function  of  0  =  a  D,  where  D  is  the  disc  major  axis.  The  results  are 
about  70%  of  those  for  a  circular  disc  of  diameter  equal  to  the  ellipse 
major  axis,  for  low  values  of  0.  As  0  increases,  this  fraction  decreases. 
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Fig.  11  -  Monte  Carlo  values  of  the  rati «  a  ...  . 
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The  radiation  absorbed  by  the  disc  edge  and  that  exiting  from  the  disc 
face  were  found  separately  in  the  elliptical  disc  calculations.  These  values 
are  shown  in  Fig.  12  as  a  function  of  (3.  We  see  that  the  greater  part  of  the 
stimulated  radiation  increase  at  high  0  goes  to  the  disc  edge,  rather  than 
to  the  disc  face.  In  high  energy  systems,  the  edge  absorber  may  thu3  have 
to  withstand  large  loadings. 

D.  Discussion 

We  have  studied  the  problem  of  fluorescence  amplification,  in  which  the 
fluorescence  emitted  by  a  material  with  optical  gain  is  amplified  by  stimulated 
emission  as  it  passes  through  the  material,  thus  increasing  the  natural 
fluorescent  loss  rate.  We  have  defined  the  ratio  of  total  loss  to  fluorescence 
M  and  the  ratio  of  stimulated  to  spontaneous  loss  A. 

Fluorescence  amplification  causes  pumping  to  become  more  and  more  difficult 
as  the  stored  energy  (and  thus  the  gain)  increases  in  a  material.  In  effect, 
the  energy  decays  with  an  instantaneous  fluorescent  lifetime  given  by  t/M, 
where  t  is  the  natural  or  low-gain  lifetime.  This  means  that  as  the  pump 
strength  is  increased,  inversion  will  not  increase  in  proportion.  The  resulting 
pump  efficiency  decrease  is  a  slow  function  of  the  pump  energy,  rather  than 
the  abrupt  upper  limit  due  to  parasitic  oscillation  (Section  III).  The  exact 
details  of  the  process  depend  on  the  geometry  of  the  material  and  the 
fluorescent  line  shape,  but  in  general  pumping  is  quite  difficult  when  the 
gain  across  the  longest  dimension  of  the  material  is  about  exp  (5)  a  150 
or  more  (we  assume  that  the  natural  lifetime  is  not  extremely  long,  so  that 
its  decrease  by  a  factor  of  10  is  a  serious  problem).  The  effect  of  fluore¬ 
scence  amplification  on  pumping  will  be  considered  in  detail  in  Section  IV. 
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12  -  Amount  of  ratio  M  of  total  to  spontaneous  loss  due 
to  radiation  absorbed  at  edge  (Mg)  and  lost  from  faces  (Mp) 
in  an  elliptical  disc,  as  a  function  of  p.  Most  of  the 
increased  loss  due  to  fluorescence  amplification  is  seen 
to  go  to  the  edge  rather  than  out  the  faces. 


The  line  shape  of  the  fluorescent  emission  (and  gain)  is  important  because 
a  line  with  broad,  low  wings  has  less  fluorescence  amplification  than  a  more 
rectangular  line,  for  the  same  peak  value  of  the  gain.  For  low  gain,  we 
have  shown  that  the  ratio  A  has  the  relative  values  2:^2: 1  for  flat  (rectangular), 
Gaussian,  and  Lorentzian  gain  profiles;  Monte  Carlo  calculations  show  that 
these  relations  are  roughly  maintained  at  higher  gains.  Actual  laser 
materials  often  have  complicated,  multi-peaked  line  profiles,  but  in  principle 
we  could  find  the  single-wave length  value  A^  for  the  geometry  we  are  using 
and  then  (perhaps  numerically)  integrate  this  with  the  amplitude  weighting 
function  corresponding  to  the  laser  material  we  are  using.  In  practice,  it  is 
undoubtedly  sufficient  to  estimate  where  in  the  rlat-Gaussian-Lorentzian 
family  the  actual  line  shape  lies  and  to  interpolate  to  find  A  for  that 
line  shape. 

The  geometry  of  the  amplifying  material  strongly  influences  the  amount  of 

fluorescence  amplification.  In  the  simple  case  of  a  sphere,  A  and  M  depend 

only  on  the  product  of  the  gain  coefficient  a  and  the  disc  diameter  D,  and 

of  course  on  the  line  shape.  When  we  consider  a  circular  disc,  we  find  that 

A  and  M  depend  not  only  on  aP,  but  also  on  the  refractive  index  ratio  between 

the  inside  and  outside  of  the  disc,  and  on  the  ratio  of  thickness  to  diameter 

of  the  disc.  An  index  match  between  inside  and  outside  is  desirable  since  it 

leads  to  the  least  possible  fluorescence  amplification.  However,  in  short- 

pulse  high-power  lasers  the  index-matching  material  may  cause  problems  due  to 

self- focusing  damage,  so  this  remedy  is  not  always  available.  The  fluorescence 

amplification  is  also  least  for  vanishingly  thin  discs,  but  of  course  such 

discs  have  insufficient  gain  along  the  path  of  the  actual  laser  pulse  to  even 
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overcome  their  own  surface  losses.  Since  the  increase  of  fluorescence 
amplification  with  thickness  is  slow,  it  is  not  harmful  to  increase  the  disc 
thickness  until  it  is  limited  by  other  considerations.  For  example,  the 
deposition  of  pump  energy  changes  with  disc  thickness,  but  this  question  is 
beyond  the  scope  of  this  paper.  An  elliptical  disc  has  the  same  type  of  index 
and  thickness  behavior  as  a  circular  disc,  but  the  fluorescence  amplification 
is  somewhat  less  than  in  the  circular  case.  Elliptical  discs  have,  the  addition¬ 
al  advantages  of  smaller  size  and  hence  lower  cost  (when  compared  to  circular 
discs  tilted  at  Brewster's  angle  to  a  circular  beam)  and  closer  coupling  to 
the  pump  sources. 

III.  PARASITIC  OSCILLATION 

Parasitic  oscillation  takes  place  when  a  material  with  optical  gain  has 
in  it  (or  through  it)  a  light  path  which  returns  on  itself.  Under  this 
condition,  the  material  will  break  into  oscillation  when  the  gain  is  large 
enough  to  overcome  the  path  losses.  Such  undesired  oscillation  will  make 
it  impossible  to  increase  the  stored  energy  in  the  oscillating  volume,  since 
once  oscillation  begins  all  further  pumping  is  immediately  converted  into 
oscillating  power  and  lost.  If  the  mode  of  oscillation  fills  an  appreciable 
fraction  of  the  laser  volume,  it  thus  sets  a  sharp  upper  limit  to  the  available 
stored  energy,  and  therefore  the  gain.  If  the  mode  gain  does  not  vary  greatly 
with  frequency  across  the  gain  line  profile  (or  if  it  oscillates  so  rapidly 
as  to  effectively  sample  all  values  of  gain  in  the  line),  then  the  place 
where  oscillation  will  begin  is  the  line  peak.  This  means  that  the  maximum 
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gain  coefficient  a  determines  the  oscillation  threshold,  so  that  rectangular 
line  profiles  are  the  best  from  the  view  point  of  parasitic  oscillation, 
since  they  maximize  the  available  laser  output  at  the  onset  of  oscillation. 

Parasitic  oscillation  cannot  be  calculated  by  ordinary  Monte  Carlo 
methods,  since  the  power  in  oscillating  modes  typically  flows  in  precise 
directions  at  any  point  in  space.  Unless  the  Monte  Carlo  algorithm  is 
lucky  enough  to  start  a  ray  at  exactly  the  correct  position  going  in  exactly 
the  correct  direction,  it  will  not  find  the  mode  at  all.  Imagine,  for 
example,  the  simple  case  of  two  reflecting  plane  mirrors  1  cm  in  diameter 
spaced  10  cm  apart,  with  a  laser  material  between  them.  In  the  case  of 
such  large  mirrors,  we  can  ignore  diffraction  and  do  a  simple  ray-tracing 
analysis.  The  mode  occupies  the  whole  volume  between  the  mirrors,  so  the 
Monte  Carlo  method  just  has  to  get  the  angle  right.  In  fact,  however,  there 
is  only  one  exactly  correct  direction  -  the  direction  perpendicular  to  the 
mirrors.  Rays  at  a  small  angle  from  the  correct  direction  will  bounce  between 
the  mirrors  a  number  of  times  before  they  miss  a  mirror  and  exit  from  the 
system.  On  the  average,  a  ray  at  an  angle  €  from  the  system  axis  will 
bounce  N  times,  where  N  is  given  by 


where  d  is  the  mirror  diameter  and  l  is  the  inter-mirror  distance.  If  the 
Monte  Carlo  method  starts  rays  at  random  in  angle,  the  probability  of  getting 
N  or  more  bounces  from  a  ray  is 
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P<N)  -  (m)  • 

The  rapid  decrease  of  p  as  N  increases  means  that  a  very  large  number  of  rays 
must  be  started  to  get  a  ray  with  many  bounces.  For  our  1  cm  by  10  cm 
resonator,  we  must  start  16,000,000  rays  (on  the  average)  to  get  one  with 
100  or  more  bounces.  Clearly,  random- starting  ray-tracing  methods  are  not 
suitable  for  mode  finding  unless  considerable  modification  is  made  to  them. 
Without  such  modification,  a  Monte  Carlo  algorithm  may  give  what  seem  to  be 
perfectly  reasonable  results  even  in  the  presence  of  an  oscillating  mode 
above  threshold. 

Since  our  Monte  Carlo  program  is  unsuitable  for  the  analysis  of 
parasitic  oscillation,  a  simplified  analytic  method  was  used  instead. 

The  method  uses  a  simple  ray- tracing  viewpoint  which  ignores  diffraction 
effects.  In  addition,  phase  is  ignored  and  only  amplitude  is  considered. 

Both  these  simplifications  are  justified  on  the  basis  that  typical  discs  used 
in  large  lasers  are  very  much  larger  than  a  wavelength,  so  that  modes  have  low 
diffraction  loss  and  close  spacing  in  wavelength. 

We  will  first  consider  lossless  modes  which  rely  on  total  internal 
reflection,  then  modes  which  have  some  loss  and  thus  require  gain  in  the 
laser  material  to  achieve  threshold,  and  finally  the  effects  of  a  rough 
disc  edge  on  the  mode  threshold. 

A.  Lossless  Modes  in  a  Circular  Disc 

It  is  possible  for  a  mode  which  is  totally  lossless  to  exist  in  a 
circular  cylinder  because  all  reflections  are  by  total  internal  reflection. 
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Consider,  for  example,  the  disc  shown  in  Fig.  13.  The  disc  faces  are  in 
contact  with  a  medium  with  refractive  index  n^,  the  disc  itself  has  index 
n^,  and  the  (smooth)  edge  of  the  disc  is  clad  with  an  absorbing  material 
of  index  (we  Ignore  any  imaginary  component  of  the  absorber  index) . 

If  the  edge  coating  index  n^  is  less  than  n^,  then  total  internal  reflection 
is  possible  off  the  inside  of  the  disc  edge,  and  a  ray  can  bounce  losslessly 
around  the  disc  perimeter  in  a  plane  parallel  to  the  faces  (Fig.  14).  In 
addition,  if  the  outside  index  n^  is  less  than  the  disc  index  n^,  a  ray 
can  bounce  from  face  to  face  of  the  disc  by  total  Internal  reflection 
(Fig.  15). 

In  general,  a  ray  may  bounce  off  both  the  faces  and  the  ed[;e,  and  the 

ray  path  will  be  complicated  (Fig.  16).  We  can  describe  such  a  path  in  terms 

of  two  angles:  Y  measures  the  angle  of  the  ray  from  a  disc  diameter  (Fig.  14), 

and  0  measures  the  angle  from  a  normal  to  the  faces  (Fig.  15).  The  angle  of 

incidence  on  the  faces  is  thus  6,  and  the  angle  0  of  Incidence  on  the  edge 

obeys  C0S<t>  ■  COS Y  sin  8.  Thus  as  0  becomes  smaller,  0  becomes  larger.  For  a 

ray  to  be  lossless,  we  must  have  both  0  >  0  =  sin  ^(n,/n„)  and 

c  12 

0  >  0c  =  sin  ^"(n^/r^).  This  dual  condition  sets  a  lower  bound  on  Y,  since 

we  must  have  COS Y  <  COSQ  /sin  0  .  Thus  we  must  have 

c  c 


which  implies  that  the  lossless  oscillation  exists  only  outside  a  center 
portion  of  the  disc  of  radius  R  ■  D/2  sin  Y  ,  where  Y  is  the  minimum  value  of  Y. 
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Fig.  13  -  The  disc  configuration  considered  in  the  parasitic 
oscillation  calculation.  A  circular  disc  of  index  n«  is 
embedded  in  a  material  of  index  m .  The  edge  of  thedisc  is 
clad  with  an  absorbing  material  of  index  n3. 


Fig.  14  -  Typical  lossless  ray  path  possible  if  the  disc  index  is 
greater  than  the  edge-coating  index.  Total  Internal  reflection 
takes  place  at  the  reflection  points. 
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Fig.  16  -  The  path  of  a  ray  which  bounces  off  both  the  faces 
the  edge.  Between  edge  reflections,  the  path  zig-zags  in  a 
plane  perpendicular  to  the  faces;  the  path  changes  to  a  new 
plane  at  each  edge  reflection. 


We  conclude  that  if  the  edge  coating  refractive  index  is  less  than  the 
disc  index  lossless  modes  can  exist  in  a  circular  disc.  This  means  that 
as  soon  as  there  is  any  gain  at  all  in  the  material  oscillation  will  commence 
and  further  gain  rise  will  be  inhibited.  The  actual  mode  volume  depends  on 
the  edge/disc  and  outside/disc  refractive  index  ratios.  It  is  clearly 
desirable  to  have  the  coating  index  larger  than  the  disc  index,  since  then 
lossless  modes  do  not  exist. 

Lossy  Modes  in  a  Circular  Disc 

Even  if  there  are  no  lossless  modes  in  a  disc,  it  will  still  oscillate  on 
lossy  paths  if  the  gain  is  high  enough.  We  wish  to  find  the  maximum  gain 

achievable  without  parasitic  oscillation,  so  we  must  find  the  mode  with  the 
least  loss. 

Consider  a  path  which  reflects  off  the  disc  faces  by  total  internal 
reflection,  but  which  has  only  partial  reflection  at  the  edge.  We  will  use 
the  same  angles  Y,  9  as  we  used  in  the  previous  section  to  describe  the  ray. 

Let  us  assume  that  the  edge  reflection  is  given  by  R(j,  the  larger  of  Fresnel's 
reflection  coefficients.  Then 


R 


sin(0Yr)J 


where  r  -  sin-1  (n^  sin  0)  and  0,  as  before,  is  given  by  0  =  m’h n^n 
OTY)  when  9  is  equal  to  0^  The  path  length  between  edge  bounces  is 
D  casV  sin  9C  -  D  n2/ni  AWY,  where  D  is  the  disc  diameter.  To  oscillate,  the 
net  gain  must  be  equal  to  unity,  or 
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e 


sin(<fl-r)  2 
sin(0hr)_ 


a  Dn2 


easy 


1 


where  a  is  the  gain  coefficient  at  the  line  peak,  0  =  «tf  1(n^/n2  cos Y)  and 
r  -  sin"^(n2/n2  sin  0).  We  could  now  vary  Y  in  this  expression  to  find  the 
minimum  a  required  for  oscillation,  but  it  is  sufficient  to  examine  the 
behavior  at  Y  ■  0  and  Y  ■  rt/2.  Solving  for  aP,  we  have 


and 


P(y)  =  a  d  = 


2n ^B/n  ■ 


n„ 


sin(0>-r) 

sin(A-r) 

easy 


\ 


where  p  -  n^/n2  and  o  *  while 


4o 

P(n/2)  a  J  2 

We  have  again  used  the  variable  p  *  aP. 

The  calculation  is  only  valid  for  n^  <  n2  and  n^  >  n2,  or  p  <  1  and 
a  >  1.  In  this  region  p(0)  is  always  less  than  P(n/2),  so  we  find  that  the 
lowest  threshold  modes  are  those  in  which  the  ray  proceeds  across  a  diameter 
of  the  disc,  bouncing  from  face  to  face  at  the  maximum  angle  which  will 
produce  total  internal  reflection.  Such  modes  obviously  fill  the  entire 
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disc  volume.  The  threshold  condition 


may  be  written 


hCRIT 


+ »,  \ 

n^v+n/.„;j  • 


Since  the  »or.t-case  ray,  are  those  that  travel  diametrically  .cross  a 

circular  disc,  we  expect  the  threshold  for  an  elliptical  disc  of  major  axis 

equal  to  the  disc  diameter  to  he  essentially  the  same  as  the  circular  threshold, 

since  the  same  path  Is  available.  The  surface  curvature  at  the  reflection 

points  Is  of  courae  larger  In  the  elliptical  case,  but  this  should  raise  the 

threshold  only  slightly,  since  all  it  does  Is  make  the  resonator  somewhat 
more  unstable. 

In  addition  to  the  mode.  Involving  reflection  off  the  edges,  a  relatively 
thick  disc  may  also  oscillate  between  It.  faces.  The  threshold  for  this  process 
is  easily  calculated,  ,i„ce  the  face  reflection  at  normal  Incidence  1, 

2 


-fo-'-M' 
\  %  +  ) 


ana  the  gain  from  face  to  face  is  lust  pxd  (ni  .  T  .  . 

just  exp  (CtL,,  e  L  is  the  disc  thickness. 

We  must  have  R  exp(ah)  =  1,  or 


"oL*CRIT  ~  2 


Of  course,  this  result  also  applies  to  elliptical  discs. 
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c •  Combined  Oscillation  Diagram 

Me  may  shov  the  effects  of  both  lossless  and  lossy  modes  on  a  single 
diagram  In  the  p,  o  plane  (recall  p  -  n,/n2  and  0-  nj/n^.  when  lossless 
modes  fill  the  disc  from  a  radius  R  outward,  we  use 

W  =  p  V1  ■  o2 

and 


to  derive 

[‘ -  (“)  ]  e2  +  o2  -  1 

Thus  the  locus  of  constant  R  in  the  p,  c  plane  is  an  ellipse.  This  result 
applies  if  p  <  1  and  o  <  1;  if  p  >  1  then  there  is  no  total  internal  reflection 
off  the  faces  and  the  answer  becomes  independent  of  p,  since  the  oscillations 

lie  in  planes  parallel  to  the  disc  faces  and  ignore  the  condition  at  the 
faces. 

Lossy  modes  obey 


so  we  may  draw  loci  of  constant  6  In  the  p,  o  plane  to  show  how  much  gain 
may  be  achieved  before  oscillation  commences.  The  lossy  mode  calculation  la 
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valid  if  o  >  1  and  p  >  1,  and  again  if  p  >  1  Che  expression  becomes  independent 
of  p. 

The  diagram  showing  lines  of  constant  oscillation- free  radius  R  and  lines 
of  constant  across-disc  gain  log  (3  Is  shown  in  Fig.  17.  From  the  view  point 
of  oscillation  suppression,  the  best  place  to  operate  is  with  both  edge  and 
face  index-matched,  since  oscillation  is  impossible  under  these  circumstances. 
If  other  considerations  make  such  operation  undesirable,  the  degree  of 
oscillation  difficulty  encountered  with  other  choices  of  refractive  indices 
is  easily  evaluated  from  the  diagram. 

D.  Parasitic  Oscillation  with  a  Rough  Edge 

Since  the  oscillatory  modes  considered  above  rely  on  rays  which  follow 
specific  paths,  one  method  of  reducing  oscillation  tendency  would  seem  to 
be  the  roughening  of  the  disc  edge,  so  that  an  initially  parallel  beam  would 
be  smeared  out  in  angle.  This  ought  to  remove  much  of  the  energy  from  the 
mode,  thus  increasing  the  loss.  However,  energy  will  also  be  scattered  into 
the  mode  from  rays  initially  not  involved.  Clearly,  a  quantitative  analysis 
of  this  situation  is  required.  Again  a  Monte  Carlo  approach  is  inapplicable, 
since  a  ray  which  is  randomly  scattered  at  the  edge  has  the  same  low 
probability  of  getting  into  the  high-gain  region  of  the  mode  as  was  the 
case  with  specular  modes.  Of  course,  we  expect  the  active  mode  volume  in 
space  and  angle  to  be  larger  with  a  rough  edge  than  in  the  specular  case, 
but  the  amount  by  which  this  helps  the  Monte  Carlo  method  is  hard  to  estimate. 

The  rough-edge  oscillation  was  therefore  analyzed  by  a  method  related 
to  the  well-known  method  of  analyzing  modes  in  an  optical  resonator.  An 


F-56 


OSCILLATION 
THRESHOLD 
GAIN  LOO 


RADIUS 
fRACTiON 
FREE  OF 
LOSSLESS 
OSCILLATION 


1.6  J- 


LOSSLESS 
OSCILLATIONS 
FILL  ENTIRE 
DISC 


Fig.  17  -  Combined  diagram  of  lossless  and  lossy  parasitic 
oscillation  as  a  function  of  the  index  ratios  of  a  circular 
disc.  In  the  lower  part  of  the  diagram  (03/02  <  1),  the 
fraction  of  the  disc  free  of  lossless  oscillations  is  shown. 
In  the  upper  part  (03/02  >  1),  the  ?  required  to  commence 
oscillation  along  lossy  paths  is  shown.  See  text. 
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initial  angular  distribution  of  optical  power  was  assumed  at  the  disc  edge, 
and  the  power  was  then  transported  across  the  disc  (with  gain),  reflected 
with  loss,  and  scattered  in  angle.  In  general,  the  distribution  after  this 
process  was  different,  in  both  amplitude  and  shape,  from  the  original 
distribution.  The  gain  coefficient  was  adjusted  after  each  iteration  to 
keep  the  amplitude  constant,  and  the  process  was  reiterated  until  the 
angular  distribution  of  power  had  reached  an  equilibrium  state.  After 
many  such  iterations,  the  gain  had  the  value  appropriate  to  the  oscillation 
threshold,  and  the  angular  power  distribution  was  that  of  the  oscillating 
mode.  Once  again,  diffraction  and  phase  were  ignored. 

Since  radiation  which  is  totally  internally  reflected  from  the  disc 
faces  has  much  higher  gain  travelling  across  the  disc  than  radiation  which 
is  only  partially  reflected,  the  power  was  assumed  totally  lost  when  it 
was  partially  reflected.  With  this  assumption,  all  points  on  the  disc 
edge  are  equivalent,  and  only  the  angle  of  the  optical  power  is  important 
(as  opposed  to  its  position).  The  hemisphere  of  possible  angles  was 
divided  into  sectors  of  equal  solid  angle,  and  the  resulting  matrix  of 
powers  was  iterated  as  above.  The  edge  roughening  was  simulated  by  taking 
the  unit  vector  along  the  direction  of  the  reflected  ray,  and  adding  to  it  a 
random  vector  with  uniform  distribution  throughout  a  sphere  of  radius 
F  COSQ,  where  F  s  1  and  9  is  the  angle  from  the  reflected  ray  to  the  normal. 
The  resulting  vector  was  then  renormalized  to  unity.  Various  amounts  of 
surface  roughening  were  simulated  by  varying  F  from  zero  (specular  reflection) 
to  one  (very  rough).  The  edge  reflectivity  R  was  independent  of  angle. 
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The  increase  in  threshold  due  to  increasing  F  is  shown  in  Fig.  18 
for  various  values  of  the  edge  reflectivity.  There  is  obviously  improvement 
over  the  specular  case,  but  it  is  not  as  great  as  might  be  expected.  From 
a  practical  viewpoint,  it  is  very  hard  to  reduce  the  reflectivity  of  a  rough 
surface  to  the  low  values  achievable  by  index-matching  a  smooth  surface, 
because  of  the  many  pits  and  cracks  in  a  rough  edge.  Therefore  one  should 
be  cautious  when  attempting  laser  improvements  by  edge  roughening. 

E.  Discussion 

We  have  shown  that  lossless  modes  exiat  for  certain  refractive  index 
ratios  between  the  inside,  outside,  and  edge  of  a  circular  disc.  Even  if 
these  modes  are  suppressed  by  proper  choica  of  the  indices,  lossy  modes 
will  still  exist.  Such  modes  place  an  upper  limit  on  gain,  since  once  the 
gain  reaches  the  mode  threshold  for  oscillation,  no  further  gain  increase  is 
possible  in  the  mode  volume.  Highest  gain  is  achieved  for  close  index 
matching  between  the  disc  and  the  material  in  contact  with  its  faces,  and 
between  the  disc  and  the  absorbing  coating  on  its  edge.  We  have  shown 
what  gain  is  possible  if  proper  matching  is  not  feasible  for  reasons  not 
connected  with  oscillation. 

Even  if  the  absorbing  material  on  the  edge  is  a  good  index  match  to  the 
disc,  there  may  be  areas  of  less-than-perfect  adhesion,  or  chips  and  cracks 
in  the  disc  edge.  Such  flaws  will  have  reflection  larger  than  the  coating- 
disc  interface,  and  in  practice  will  make  the  achievement  of  very  high 
absorption  at  the  edge  difficult.  If  the  higher-reflectivity  flaws  are  small, 
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diffraction  will  spread  the  reflected  energy,  but  we  have  seen  that  such 
spreading  causes  only  small  threshold  improvements.  Thus  a  well-index- 
matched  material  must  be  applied  extremely  well  to  a  flawless  surface  if 
a  high  parasitic  oscillation  threshold  i*  to  be  achieved. 

IV.  EFFECTS  ON  PUMPING 

Our  discussion  so  far  has  been  confined  to  an  instantaneous  snapshot 
of  the  fluorescence  amplification  and  parasitic  oscillation  processes. 
Lasers  are  usually  pumped  by  relatively  slow  pump  pulses,  and  we  are 
interested  in  knowing  what  decrease  in  efficiency  is  caused  by  the  loss 
processes  we  have  calculated  during  such  a  pulse.  We  will  first  find  the 
effect  of  fluorescence  amplification,  and  then  the  effect  of  parasitic 
oscillation. 

A.  Fluorescence  Amplification 

In  the  low-gain  case,  when  fluorescence  amplification  is  negligible, 
the  stored  energy  decays  with  -the  fluorescent  lifetime  t.  The  gain  has  the 
same  behavior,  and  so  we  may  write  a  differential  equation  for  the  across- 
disc  gain  log  p  in  the  form 


P  -  p(t)  -  0/t 

where  p(t)  is  the  pump  pulse  (in  appropriate  units).  At  higher  gain, 
fluorescence  amplification  increases  the  loss  by  a  factor  M,  and  so  we  have 

P  -  P(0  -  P  M  (P)/T. 
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If  we  normalize  to  t  •  1,  and  assume  a  half-sine  approximation  to  the 
pulsa,  we  get 


pump 


•  Jtp  /  V 

P  “  ~2T  8ln  “  P*CP)  0  «  t  «  T 

where  pQ  i8  the  area  of  the  pump  pulse  and  T  is  the  base  width  of  the  pump 
pulse  (in  units  of  the  fluorescent  lifetime).  In  the  range  0  s  M  £  5,  we 

may  make  a  good  approximation  to  M  for  any  of  the  geometries  and  lineshape. 
studied  by  using  the  form 


M  s-  exp  p  +  G  p2  +  H  p3j  . 

Note  that  F  may  be  found  from  the  low-gain  value  of  A,  since  A  &  F  £  when 

0  «  1.  We  may  then  find  G  and  H  by  fitting  a  straight  line  to 
(fin  M/p  -  F)/p  «*.  G  +  Hg. 

Once  M  is  approximated,  we  may  integrate  the  differential  equation  for 
3  until  it  reaches  its  peak  value.  A  plot  of  such  maximum  values  is  shown 
in  Fig.  19  as  a  function  of  the  pumping  pQ,  for  various  values  of  the  pump 
pulse  width.  We  see  that  fluorescence  amplification  is  more  serious  for 
Bla">  lar8e"ar«a  P»®P  pulses  while  faster  or  smaller  pulses  suffer  less 
loss.  In  the  design  of  a  flashlamp-pumped  laser,  the  effects  of  current 
density  on  the  pumping  efficiency  of  flashlamps,  and  the  effects  of  energy 
input  on  flash  lamp  life,  should  be  integrated  with  this  decay  information 
ia  order  to  choose  the  optimum  operating  point. 

Spatially  uniform  gain  has  been  assumed  in  all  our  calculations  of 
fluorescence  amplification.  However,  the  removal  of  energy  due  to 
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fluorescence  amplification  will  be  spatially  non-uniform  in  general.  Thus  the 
gain  will  become  non-uniform  during  a  pumping  pulse  even  if  the  pumping  is 
uniform.  The  analysis  of  this  problem  is  beyond  the  scope  of  this  paper,  but 
it  is  clear  that  proper  tailoring  of  the  spatial  distribution  of  the  pump 
energy  can  reduce  or  eliminate  the  difficulty. 

B.  Parasitic  Oscillation 

Parasitic  oscillation  causes  a  rapid  limitation  of  stored  energy,  rather 
than  the  slow  limitation  we  have  seen  above  for  fluorescence  amplification. 
This  is  because  a  mode  below  its  oscillation  threshold  contributes  little 
to  the  loss  rate,  but  once  the  mode  is  over  threshold  the  loss  rate  is  equal 
to  the  energy  input  rate.  Thus  the  gain  is  sharply  limited  at  the  oscillation 
threshold.  Figure  20  shows  again  the  curve  of  gain  versus  pump  area  for  a 
pump  pulse  of  base  width  equal  to  the  fluorescent  lifetime,  and  also  shows 
the  type  of  gain  saturation  expected  from  the  onset  of  oscillation.  The 
oscillation  threshold  is  shown  for  various  ratios  of  the  refractive  index 
of  the  disc  edge  coating  to  the  refractive  index  of  the  disc;  this  ratio 
determines  the  reflectivity  at  the  edge.  We  have  again  assumed  a  disc-to- 
environment  index  ratio  of  1.56.  This  diagram  illustrates  that  unless 
very  good  index  matching  (or  other  reflection  suppression)  is  used,  parasitic 
oscillation  will  limit  available  energy  storage  well  before  fluorescence 
amplification  losses  become  appreciable. 

V.  CONCLUSIONS 

The  effects  of  fluorescence  amplification  and  parasitic  oscillation  on 
laser  energy  storage  and  pumping  efficiency  have  been  analyzed.  Fluorescence 
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amplification  .at,  no  definite  npper  11. It  on  stored  energy,  but  Instead 
makes  punping  note  and  note  difficult  as  the  energy  density  rises.  Parasitic 

oscillation,  on  the  other  hand,  set.  an  ahrupt  upper  Unit  to  stored  energy 
at  the  oscillation  threshold. 

fluorescence  anpllflcation  nlninlced  h,  Index-matching  the  naterlal 
in  contact  with  the  faces  of  a  laser  disc  to  the  refractive  Index  of  the 
hlec.  increasing  the  disc  thickness  cause,  only  .  slow  lnctease  In  the  anount 
of  fluorescence  anpllf-. cation,  .0  thicker  disc,  are  of  value  if  p„np  unlfornlty 
aw?  thermal  distortion  limits  are  not  exceeded. 

Parasitic  oscillation  will  occur  at  very  low  gain  level,  if  the  edge 
cladding  naterlal  ha,  a  lower  refractive  Index  than  the  disc  naterlal.  If 
this  error  is  avoided,  the  level  at  which  oscillation  occur,  will  he 
maximized  If  the  cladding  Index  Is  close  as  possible  to  the  disc  Index. 

In  addition,  Index-natchlng  the  faces  raise,  the  oscillation  threshold. 
Roughening  of  the  edge.  doe.  not  cause  a  great  increase  of  the  oscillation 
threshold,  and  nay  increase  the  average  reflectivity  of  the  edge,  thus 
doing  more  harm  than  good. 

Khat  imitation,  do  fluorescence  amplification  and  parasitic  oscillation 
Place  on  glass  laser  amplifier  disc  size?  Figure  20  show,  that  parasitic 
oscillation  1,  the  worst  prohlen.  Unless  stringent  precaution,  are  taken, 
oscillation  will  set  In  ,t  an  acr.ss-dlsc  gain  of  exp(3)  „t  less,  assuming  th. 

dUC  l”  “  *VOld  3elf" focusing  problems,  with  .5  J  cm-3  stored,  the 
Une-peak  gain  will  be  on  the  order  of  .1  Thl,  Implies  .  3„  „  disc 

•t  most.  This  value  will  vary  .ith  stored  energy  and  the  „»unt  of  gain  per 
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unit  stored  energy,  but  values  much  smaller  than  those  quoted  increase 
the  amount  of  energy  lost  to  the  fixed  loss  coefficient  in  the  glass. 

Thus  operation  with  discs  of  dimension  larger  than  about  30  cm  will  be 
very  difficult,  and  operation  in  the  30  cm  region  will  require  great  care 
in  edge  treatment  to  suppress  parasitic  oscillations. 
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ABSTRACT 


EY-1  is  a  terbium-silicate  glass  recently  developed  by 
Owens-Illinois.  It  offers  great  potential  as  a  Faraday  rotator- 
isolator  in  high  power  laser  systems  because  of  its  high  Verdet 
constant,  V  =  -4.06  ±  0.1  x  10"2  min/Oe-cm  at  1.064  m«  This  allows 
one  to  use  shorter  lengths  of  glass  and  lower  magnetic  fields  than 
with  other  rotator  glasses.  We  have  measured  V  in  this  glass  at 
1.06  m  with  fields  up  to  60  kOe  and,  in  addition,  the  change  in  V 
with  temperature  and  wavelength: 


dV\  -4 

«kT J  -  "2,2  ±  0.4  x  10  min/0e-cm-°C 

0  '1.064m 


=  1.1  ±  0.2  x  10  5  min/ Oe-cm-cm 

The  requirements  that  these  numbers  impose  on  the  field,  temperature 
and  wavelength  stability  and  uniformity  are  discussed. 


PROBLEM  STATUS 


This  is  an  interim  report  on  a  continuing  problem. 
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NRL  Problem  K03-08.502 
ARPA  Order  2062 


G-3 


I.  INTRODUCTION 


Faraday  rotators  are  used  in  laser  plasma  experiments  to  isolate 
the  laser  from  light  back-reflected  from  the  target.  These  back- 
reflections  are  especially  serious  in  large  multistage  laser  systems 
because  low  reflected  light  intensities  can  be  amplified  to  destruc¬ 
tive  levels  by  propagating  backwards  through  the  amplifying  chain. 

The  use  of  Faraday  rotators  and  polarizing  optical  elements  offers 
one  method  of  preventing  this  destructive  feedback.  This  paper 
considers  some  design  concepts  concerning  Faraday  rotators  for  laser 
systems  and  presents  measurements  made  on  Owens-Illinois  EY-1 
Faraday  rotator  glass. 

II.  DESIGN  CONCEPTS 

, - . -i?1®  design  of  a  Faraday  rotator  system  for  a  high  power  laser 
(>10  W)  is  governed  by  the  diameter  of  the  laser  beam  and  the 
length  of  the  Faraday  rotator  material.  The  length  is  determined  by 
the  self- focusing  length  of  the  material  at  the  desired  laser  power 
density.  The  Verdet  constant  and  the  length  determine  the  necessary 
magnetic  field  strength.  J 

Glasses  are  presently  used  as  the  Faraday  rotator  material  in 
high  power  laser  systems.  Recently  Owens-Illinois  has  developed  for 
commercial  sale  a  terbium- silicate  base  glass,  EY-1,  that  has  a 
Verdet  constant  of  about  -4  x  10"2  min/Oe-cm.  It  is  a  paramagnetic 
glass.  A  similar  glass  composition  was  first  formulated  for  Faraday 
rotator  purposes  by  Robinson  and  Graf.*1*  For  comparison”'3}  at 
1.06  Schott  SF-6  lead-silicate  glass  has  a  Verdet  constant  of 
+  1.9  x  10  min/Oe-cm  and  a  common  soda  lime  silicate  glass  has  a 
Verdet  constant  of  +  0.65  x  10'2  min/Oe-cm.  These  last  two  glasses 
are  both  diamagnetic.  The  high  Verdet  constant  of  EY-1  allows  a 
shorter  length  of  material  to  be  used  for  a  given  field.  This 
reduces  the  energy  in  the  magnetic  field  and  increases  the  self- 
focusing  damage  threshold. W  The  base  glass  is  similar  to  the  Owens- 
Illinois  laser  glass  ED-2  and  consequently  it  should  exhibit 
desirable  surface  and  bulk  damage  characteristics. 

III.  EY-1  GLASS  MEASUREMENTS 

We  have  measured  the  following  properties  of  EY-1  glass: 

:  constant  V  up  to  60  kG  at  1.06  u  and  19° C. 

change  in  V  with  temperature  at  constant  wavelength 
at  30  kG 

change  in  V  with  wavelength  at  constant  temperature 
at  30  kG. 
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fhese  numbers  determine  the  field  strength  of  the  system  and  the 
temperature  and  wavelength  stability  required  to  maintain  a  given 
extinction  ratio  with  crossed  polarizers. 

The  experimental  arrangement  is  shown  in  Fig.  1.  The  Verdet 
constant  was  measured  by  noting  the  change  in  rotation  with  magnetic 
field  9  -  Vtt  .  0  is  traditionally  measured  in  minutes,  t  in  cm, 

’ JT0nsequently^  ln  The  magnetic  field  was 

p  ovided  by  the  NRL  magnet  facility  and  was  constant  to  ±1%  over 

2"  VOir  dia  x  3-9  “  lon8) •  The  rotation  was  detected 

by  measuring  the  null  in  transmitted  light  through  Glan-Thompson 

prisms  with  a  photomultiplier  detector.  The  extinction  ratio  was 
100:1,  and  the  angular  error  was  ±  20  min.  The  temperature  was 
monitored  with  an  iron- cons tantan  thermocouple  that  was  thermally 
attached  to  the  glass  surface  with  silicone  grease.  The  reference 
junction  was  held  at  a  constant  temperature  in  ice  water,  and  the 

experiments^eratUre  W8S  t0  ±1#C  dUrlng  the  constant  temperature 

A  Chromatix  Model  1000C  Nd: YA-fG  laser  was  used  as  the  light 
Th®  wavelength  was  varied  by  adjusting  the  wavelength 
selector  prism  on  the  laser  to  the  different  Nd3*  laser  transitions, 
and  the  laser  was  operated  at  about  10  pps  in  a  pulsed  mode.  No 

Tdt  <°  syn^hronou8ly  d«tect  the  pulses  or  to  use  boxcar 

■nffP?tl^teCiniq‘|eB  because  the  accuracy  of  the  experiments  was 
sufficient  for  Faraday  rotator  design  purposes. 

III-l.  Verdet-Constant  vs  Field  Strength 

Figure  2  shows  the  rotation  as  function  of  field  strength. 

The  Verdet  constant  at  1.064  /<  and  19°C  is  V  -  -4.06  ±  0.1  x  10"2  min/ 
Oe-cm.  This  is  lower  than  the  value  quoted  in  Ref.  (1)  because  of 
the  lower  Tb  content  W. 

Ill- 2.  Verdet  Constant  vs  Wavelength 

.  ligU^  3rfhfVf8  the  varlatlon  ^  the  Verdet  constant  with  wave- 
■®n85h*  Tb®1!J®ld  waB  held  constant  at  30  kOe  ±  1%  and  the  temperature 
as  eld  at  12  C  t°  ±  0.2°C.  The  laser  was  tuned  over  the  available 

determine!!8  JT  i  4  t0  1,079  4’  and  th®  chanSe  in  station  was 
determined.  The  change  in  Verdet  constant  is: 


M 

^  T  =  12° C 


1.1  ±  0.2  x  10  5  min/ Oe-cm- cm  3. 
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The  rotation  increases  with  frequency  and  the  variation  per  cm“l  is 
about  10"3  x  Verdet  constant  at  1.06  /Li- 

III-3.  Verdet  Constant  vs  Temperature 

(I1-??16  4  shows  the  variation  in  the  Verdet  constant  with  tempera¬ 
ture  »  This  measurement  is  the  least  accurate  of  those  made.  The 
sample  was  heated  with  hot  air  and  allowed  to  cool  to  the  equilibrium 
magnet  temperature.  During  the  cooling  period  the  temperature  was 
found  to  be  reasonably  uniform  in  the  sample  by  changing  the  direction 
of  hot  air  flow  and  remeasuring  Also  the  rate  of  cooling  to  the 

plexiglass  holder  was  slower  than  the  calculated  rate  of  thermal  equili¬ 
brium  in  the  glass  itself  (~  5  sec) ,  which  indicates  that  a  reasonably 
flat  temperature  profile  existed  across  the  sample.  The  measurements  at 
the  highest  temperature  are  less  accurate  as  they  were  taken  under  the 
most  non-uniform  heat  distributions  in  the  sample.  The  value  is: 

|  =  -2.2  ±  0.4  x  10  4  min/0e-cm-°C. 

X  =  1.064 

This  value  per  degree  Kelvin  is  about  10”2  x  Verdet  constant  at  1.06  ju. 

The  Faraday  rotation  by  a  paramagnetic  material  as  a  function  of 
field  and  temperature  is  given  by'3*'); 


0  -  kl  tanh  (/i  H/kT)  +  BH 

SeH  (1) 

—  c  -tH/T  for  ~  «  1 

where  -  -  Verdet  constant  at  the  specified  wavelength  X,  n  is  magnetic 
moment  in  Bohr  Magnetons, H  is  field  in  Oe,  k  is  Bo ltzman's® const ant, 
l,  is  the  sample  length,  and  A  is  a  constant.  B  is  a  temperature 
independent  term  which  is  usually  found  to  be  negligible*?).  It  will 
not  be  considered  further  here.  Thus 

d  dV  _  c _  V 

dT  IT  “  2  “  “  T 

which  in  our  case  for  T  ■  12° C,  yields: 


dV 

dT 


-1.4  x  10-4 


min/0e-cm-°C. 
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This  is  less  than  the  measured  value,  and  outside  of  the  experimental 
error  limits.  A  systematic  error  may  be  due  to  the  fact  the  center 
of  che  glass  rod  is  slightly  higher  in  temperature  than  the  edge  where 
the  temperature  was  monitored.  This  is  caused  by  the  cylindrically 
symmetric  temperature  profile  of  a  cooling  thin  cylinder.  This  however 
would  not  reduce  the  value  by  more  than  20  -  30%.  We  use  the 
experimental  value  which  is  higher  than  the  theoretical  value.  The 
measured  value  allows  one  to  calculate  the  required  temperature 
uniformity  conditions  in  a  Faraday  rotator  system. 

(2  g  9) 

IV.  EXTINCTION  RATIO  VS  CHANGE  IN  VERDET  CONSTANT  *  ’ 

For  transmission  of  unpolarized  light  through  crossed  polarizers, 
which  we  assume  is  the  normal  operating  condition  of  the  Faraday 
rotator,  the  following  equation  is  obtained (see  Fig.  5): 

IT  =  k  IQ  sin2 (9  -  Jt/4)  (3) 

in  which  0  =  9(y)  where  y  can  be  field,  temperature,  or  wavelength, 
and  9  is  defined  such  that  I  =  0  at  0  *  it/4.  The  rotation  of  the 
direction  of  polarization  by  the  glass  is  nominally  set  to  be  45°, 
but  variation  in  the  rotation  with  temperature  and  wavelength  may 
occur.  To  determine  the  effect  of  these  variations  we  note  that: 

9  *  n/4  +  A0 
sin2  (0  -  :./4)  =■  A02 

A9  =  |p  Ay  (A9  measured  in  radians)  (4) 

*T  “  4  *o  (#  ^ 

Y1  =  %  Ay)  15  Extinction  Ratio  =  E.  (5) 

Using  the  relation  9  =  t  HV  where 


for 


then 


then 


19 

3y 


l  H 


2V 

ay 
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(1)  Transmission  vs  change  in  temperature  (y  ■  T  °K)  using  equation  (fi) 
and  the  measured  values  for  V  and  ■*¥,  (0  ■  45°) : 


0  =  2700 
V  =  .04 


6.74  x  10 


Oe-cm 


IV 

&T 


■2.2  x  10  4  min/0e-cm-°C  =  -6.5  x  10~®  rad/0e-cm-°C 


Y  -  0.96  x  10-5  (AT)2,  £T  in  °C. 
o 


-4 

To  maintain  an  extinction  ratio  of  <  10  ,  one  Mast  hold  the 

temperature  to  ±  1.5°C. 


(2)  Transmission  vs  change  in  wavelength  (y  - 
equation  (6)  and  the  measured  values  for  V  and 


& 


cm  1)  using 


J 

V 


6.74  x  104 


Oe-cm 


■  1.1  ±  0.2  x  10  ^  min/Oe-cm-cm  1 

-9  -1 

=  3.3  x  10  rad/Oe-cm-cm 

~  =  2.45  x  10  S  (AX)2,  AX  in  cm"1, 
o 

3+  -1  1T  -4 

In  Nd  glass,  where  AX  °*  100  cm  ,  —  °*  2  x  10  .  This  extinction 

ratio  is  tolerable  in  a  large  Nd  glass  system. 

(3)  Transmission  vs  change  in  magnetic  field  (y  -  H  in  Oe).  From 
equation  (5): 


E 


2 
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using 


as 

3H 


-  IV, 


then 


E  -  %  l2V2  AH2 


or 


Using 


H  =* 


e_ 

•tv 


4H.  M 
h  e  * 


For  an  extinction  ratio  E  of  10"^,  and  8 


45°  -  .79  rad. 


-  1.8  x  10~2  -  1.87.  . 


The  magnetic  field  inhomogeneity  must  be  less  than  1.87.. 


V.  ABSORPTION 

Faraday  rotator  glasses  typically  have  higher  intrinsic 
absorption  coefficients  than  normal  optical  glass.  This  is  due  to 
the  strong  optical  transitions  in  the  UV  and  blue  spectral  regions 
which  provide  the  field  dependent  indices  of  refraction.  A  figure  of 
merit  that  has  been  used  for  Faraday  rotator  materials^)  is  V/a. 

Where  V  is  the  Verdet  constant  in  min/Oe-cm  and  a  is  the  intrinsic 
absorption  coefficient  in  cm"1.  At  1.06  n,  a  for  EY-1  is  less  than 
0.005  cm"1,!®'  for  Schott  SF-6  it  is  w  0.001  cm'M3)  The  Schott  glass 
is  probably  superior  to  the  EY-1  glass  using  this  figure  of  merit. 
However,  the  selection  criteria  for  Faraday  rotator  glass  that  is  to 
be  used  in  high  energy  laser  systems  is  maximum  damage  resistance  and 
highest  Verdet  constant.  Residual  intrinsic  absorption  and  tempera¬ 
ture  and  wavelength  stability  are  secondary  considerations.  EY-1 
appears  to  be  a  very  suitable  material  for  high  power  optical  rotators. 
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Fig.  1  -  Experimental  arrangement 


ROTATION  THROUGH  POLARIZER  IN  DEGREES 


ROTATION  {DEO} 


Fig.  4  -  Rotation  vs  temperature  at  constant  wave¬ 
length  and  field.  H  =  30  kOe  and  X  =  1.064  u.  R>ints 
and  crosses  are  different  data  sets.  The  thermo¬ 
couple  response  is  18.9°  C/mVolt. 
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APPENDIX  H 


MODE  LOCKING  OF  CaLaSoap:Nd*t 

J.P.  Letellier 

Laser  Physics  Branch 
Optical  Sciences  Division 

Naval  Research  Laboratory 
Washington,  D.C.  20390 

ABSTRACT 


Nonlinear  effects  have  generally  precluded  the  generation  of 
time-bandwidth  limited  pulses  shorter  than  15-20  psec.  In  this  letter 
the  generation  of  5  psec  time  bandwidth  limited  pulses  from  a 
crystalline  neodymium  host  with  properties  intermediate  between 
Nd:YAG  and  Nd: glass  is  reported. 


Mode  locking  has  been  obtained  in  a  number  of  neodymium  hosts, 
but  reliable  time- bandwidth  limited  operation  has  been  obtained  only 
in  neodymium: YAG'  **' .  In  neodymium  glass,  very  early  weak  pulses  may 
exhibit  time  bandwidth  limited  behavior  and  picosecond  duration,  but 
self- focusing  and  the  consequent  self-phase  modulation  destroy  the  time- 
bandwidth  limited  character  of  the  more  energetic  pulses  in  the  pulse 
train  (as  has  been  shown  by  the  elegant  work  of  Carmen  et  al'^'). 

In  order  to  generate  time-bandwidth  limited  pulses  shorter  than 
the  15-20  psec  reported  from  Nd:YAG,  a  host  material  with  a  smaller  ^ 
value  than  Nd: glass  and  a  broader  bandwidth  than  Nd:YAG  is  required. 

This  letter  reports  the  demonstration  of  apparent  time-bandwidth 
limited  operation  with  one  such  material,  neodymium  doped  calcium 
lanthanum  silicate  oxyapatite  (CaLaSoap: Nd) .  The  properties  of  this 
material  have  been  reported  by  its  developers;  those  relevant  to  this 
investigation  are  that  the  fluorescence  linewidth  is  53  cm"^  and 
the  induced  emission  cross  section  iB  =*  1.6  x  10"^  cm^.  (4) 

A  laser  rod  5.4  cm  long  and  6.4  mm  diameter  of  CaLaSoap  doped 
with  2  wt7«  neodymium  was  flashlamp  pumped  and  mode  locked  in  a  110  cm 
long  cavity  formed  by  one  1007.  reflectivity  5  meter  radius  mirror  and 
one  flat  757.  reflectivity  mirror.  The  rod  had  flat  AR  coated  ends,  but 
had  to  be  cocked  at  two  degrees  to  the  optical  axis  to  prevent  the 
formation  of  satellite  pulse  trains.  Eastman  9740  dye  in  1,2-Dichoro- 
e thane  was  used  in  a  flowing  dye  cell  *  0.13  mm  thick. 


*  This  work  supported  by  the  Advanced  Research  Projects  Agency  under 
ARPA  Order  2062. 

t  To  be  submitted  to  Applied  Physics  Letters  for  publication. 
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T  .  ..  F1r*  1  ahorf  “  oscillo8rara  the  pulse  train  taken  from  a 
Tektronix  519  oscilloscope  at  50  nanoseconds/ cm.  The  total  energy  was 
~25  mJ  in  a  TEM00  spatial  mode.  The  pulse  train  duration  was  only 
roughly  optimized  to  be  as  short  as  possible  using  a  757.  mirror.  A 
557.  mirror  would  have  been  a  more  optimum  choice  in  principle™* .  but 
was  not  available  at  the  time.  It  should  be  noted  that  scatter  from 

the  AR  coatings  apparently  caused  some  interpulse  noise  to  develop  late 
in  the  pulse  train.  r 


The  spectrum  of  the  entire  pulse  train  was  recorded  on  a  Baird 

AJiL2’meteL8peCtr 0,116 ter  havln®  a  dispersion  figure  of  =“3.2  A/mn 
at  1.06  m.  Figure  2  shows  a  densitometer  scan  of  a  representative 
spectrogram.  The  center  frequency  was  determined  to  be  1.06113  urn 
and  the  width  at  half  intensity  was  determined  to  be  1.1  an  ±  0  2  cm"] 
No  evidence  of  gross  spectral  broadening  corresponding  to  self-phase  ' 
modulation  was  observed,  hence  the  inter pulse  noise  appeared  to 
correspond  to  scattered  light  from  the  rod  ends  with  the  same  spectral 
content  as  the  main  pulse  rather  than  to  pulse  breakup  due  to  self 
tocuslng. 

.  in  The1pulse  duration  was  examined  using  two  photon  fluorescence  in 
a  K)  cm  long  cell  of  Rhodamine  6B.  A  half  width  of  ~  5  psec  was 
obtained  by  this  method  with  a  contrast  ratio  of  -  2:1.  The  experi¬ 
mental  set  up  was  similar  to  that  of  Clobes  and  Brienza(1*  in  which  a 

W6S  obtalned  Photographically  while  a  thin  cell 
and  photomultiplier  set  up  gave  a  2.9:1  ±  .2  contrast  ratio. 

.  ,°n6/,1?al  exPerlment  was  done  in  which  the  spectra  was  found  to 
■wb.uati.lly  overlap  the  fluorescence  line  of  the  .ecord  ,tronge.t 
^transition  («f3/2  -  *Iu/2)  at  1.06111  cm  to  an  extent 
sufficient  to  permit  amplification  of  this  pulse  in  Nd: YAG  amplifiers. 
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